flD-ftl55  230 
UNCLASSIFIED 


NOTES  ON  ATMOSPHERIC  STABILITV  TURBULENCE  AND  ITS  1/1 

NEASUREMENTS<U)  ARMV  DUGUAV  PROVING  GROUND  UT 
C  A  BILTOFT  HAR  85  OPG-TR-85-204 

F/G  4/1  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

national  BURCAu  of  STANDAROS  l%i  A 


AD-A155  230 


TECOM  No.  7-C0-RD4-DP0-007 
DPG  Document  No.  DPG-TR-85-204 


TECHNICAL  REPORT 

NOTES  ON  ATMOSPHERIC  STABILITY,  TURBULENCE, 
AND  ITS  MEASUREMENT 


CHRISTOPHER  A.  BILTOFT 


Meteorological  Branch 
Test  Design  and  Analysis  Division 
Materiel  Test  Directorate 


>- 

a. 

c- 

CJ> 

Lu] 

iIZ 


US  ARMY  DUGWAY  PROVING  GROUND 
DUGWAY,  UT  84022 

MARCH  1985 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


85  5  20  19  9 


Disposition  Instructions 

Destroy  this  report  when  no  longer  needed.  Do  not  return  to  the  originator. 


Disclaimer  Statement 

The  views,  opinions,  and/or  findings  in  this  report  are  those  of  the  author(s) 
and  should  not  be  construed  as  an  official  Department  of  the  Army  position, 
unless  so  designated  by  other  official  documentation. 


Trade  Names  Statement 

The  use  of  trade  names  in  this  report  does  not  constitute  an  official  indorse¬ 
ment  or  approval  of  the  use  of  such  commercial  hardware  or  software.  This 
report  may  not  be  cited  for  purposes  of  advertisement. 


•Ificer.sion  For 

•TTS  GK/UT 
DTiC  lAB 
Unannaunofld 
Justlflcatioo. 


.4 


□ 

n 


By - 

Distribution/ 

Availability  Cod  at  j 
I  Avail  extc/oT 
|Dls^  Spaclal  I 


UNCLASSIFIED 


security  classification  of  this  page  '»Ti*n  D«(»  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCrrONS 

BEFORE  COMPLETING  FORM 

4.  TITLE  fand  Subtitle) 

NOTES  ON  ATMOSPHERIC  STABILITY,  TURBULENCE,  AND 

ITS  MEASURMENT 

5  Type  of  REPORT  4  PERIOD  COVERED 

Technical  Report 

6.  PERFORMING  ORG.  REPORT  NUMBER 

DPG-TR-85-204 

7.  AUTHOHf*; 

Christopher  A.  Biltoft 

8.  CONTRACT  OR  GRANT  NUMBERi'*; 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Dugway  Proving  Ground 

ATTN;  STEDP-MT-DA-M 

Dugway,  UT  84022 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 

3552515201 

1  1  CONTROLLING  office  NAME  AND  ADDRESS 

US  Army  Dugway  Proving  Ground 

Dugway,  UT  84022 

12.  REPORT  DATE 

March  1985 

13.  NUMBER  OF  PAGES 

56 

14  monitoring  agency  name  4  ADDRESS^//  dUierent  from  Controlling  Office) 

15.  security  class,  (of  thie  report) 

UNCLASSIFIED 

15a.  DECLASSIFI  cation/ DOWN  OR  ADI  NO 

schedule 

15  DISTRIBUTION  ST  ATEMEN  T  Cof  th/«  RepofO 

Approved  for  public  release;  distribution  unlimited. 

17  DISThiBU'^ION  statement  (of  the  abafrnct  entered  in  Block  20,  if  different  from  Report) 

10.  SUPPLEMENTARY  NOTES 

IS  KEY  WORDS  fConllnue  on  reveree  aide  if  neceatary  and  identify  by  block  number) 

Richardson  Number  friction  velocity  stability 

Obukhov  length  roughness  length 

Gaussian  diffusion  similarity 

K-theory  turbulent  kinetic  energy 

turbulence  meteorological  instrumentation 

20  A9STRACT  CCootlaua  an  rararma  ef'^  ft  naaaaaaey  aad  idarUiiy  by  block  numbat) 

Modernization  of  Dugway 's  diffusion  test  capabilities  include  the  pro¬ 
curement  of  state-of-the-art  instrumentation  and  an  upgrade  in  data  collec¬ 
tion  methodology  and  modeling.  Dispersion  of  a  buoyant,  passive  plume  is  due 
to  transport  by  mean  atmospheric  motions.  Although  K-theory  and  Gaussian 
diffusion  equations  are  burdened  by  restrictive  assumptions,  these  equations 
offer  the  most  practical  means  of  modeling  the  effects  of  turbulence  on  dif¬ 
fusing  material. 

DO  ,  1473  EWTIOt*  OF  t  HOW  6S  IS  OBSOLETE  UNCLASSIFIED 


security  classification  of  This  PA-.E  fWwn  Pete  Entered) 


UNCLASSIFIED 

security  classification  of  This  PAGgfWi«n  Dmia  Entmnd) 

BLOCK  20.  ABSTRACT  (cont'd) 

Turbulence  occurs  as  vortex  stretching  in  a  sheared  environment  causes 
an  energy  cascade  towards  progressively  smaller  eddy  sizes  until  viscous 
dissipation  occurs.  Direct  turbulence  measurements  require  an  assortment  of 
instrumentation  due  to  the  range  of  eddy  sizes.  Consequently,  it  is  most 
economical  to  estimate  turbulence  from  stability  parameters  and  most  models 
accept  such  turbulence  estimators  as  input.  The  Obukhov  Length,  a  stability 
parameter  obtained  from  principal  turbulent  kinetic  energy  components,  can  be 
obtained  from  tower  profile  data  using  flux-gradient  relationships. 

LED  scintillometers,  acoustic  sounders,  radiometers,  precision  thermo¬ 
meters,  and  sonic  anemometers  are  needed  to  upgrade  Dugway's  turbulence  mea¬ 
surement  capabilities.  Data  from  these  instruments  will  help  define  test  go/ 
no-go  conditions,  improve  quality  control,  enhance  posttest  analysis  capabil¬ 
ities,  and  provide  improved  input  for  diffusion  modeling. 


_ UNCLASSIFIED _ 

security  classification  of  this  PAGEriFhen  Data  EnlaradJ 


NOTES  ON  ATIOSPHERIC  STABILITY,  TURBULENCE,  AND  ITS  MEASUREMENT 

TABLE  OF  CONTENTS 

ACKNOWLEDGEMENTS . 

EXECUTIVE  SUMMARY . 

1.  INTRODUCTION . 

2.  ^  DEFINITION  OF  TURBULENCE . 

3.  THE  LOGARITHMIC  AND  POWER  LAW  WIND  PROFILES... 

4.  K-THEORY  AND  GAUSSIAN  DIFFUSION . 

5.  KINETIC  ENERGY  AND  THE  RICHARDSON  NUMBERS . 

6.  THE  OBUKHOV  LENGTH  AND  z/L . 

7.  STABILITY,  TURBULENCE,  AND  DIFFUSION . 

8.  MEASUREMENT  REQUIREMENTS  AND  INSTRUMENTATION.. 

9.  CONCLUSIONS . 

APPENDIX  A.  SYMBOLS . 

APPENDIX  B.  REFERENCES . 

APPENDIX  C.  ABBREVIATIONS . 

APPENDIX  D.  DISTRIBUTION  LIST . 


ACKNOWLEDGEMENTS 


I  wish  to  express  my  appreciation  for  the  thorough  technical  review  of 
this  document  by  Mr.  John  White  of  the  Dugway  Meteorology  Branch  and  by  Mr. 
James  and  Ms.  Christina  Wheeler  of  the  Dugway  Technical  Editing  Section.  My 
thanks  are  also  extended  to  Ms.  Shauna  Edgeman  of  the  Word  Processing  Section 
for  performing  nontechnical  editing  and  text  revision. 


I 


I 


I 


EXECUTIVE  SUMMARY 


This  :’)cument  examines  some  basic  relationships  between  atmospheric  dif¬ 
fusion,  tur  )ulence,  and  stability.  It  lists  the  measurements  required  for  an 
adequate  de  .cription  of  the  atmosphere  for  diffusion  applications  and  the 
types  of  eq'iipment  needed  for  obtaining  these  measurements. 

A  substance  released  into  the  atmosphere  is  dispersed  by  atmospheric 
motions.  Turbulent  motions  determine  the  rate  of  spread  of  the  material. 

Mean  motions  transport  the  material  downwind.  Dispersion  involves  scales  of 
motion  ranging  from  the  wavelengths  of  light  (measured  in  micrometers)  to 
synoptic-scale  motions  (measured  in  many  kilometers). 

Creation  of  turbulence  begins  when  the  incoming  flux  of  solar  radiation 
is  absorbed  by  the  earth's  surface  and  part  of  it  is  reradiated  as  heat.  This 
causes  thermal  discontinuities  near  the  ground.  These  discontinuities  and 
friction  from  the  interaction  of  surface  roughness  elements  and  viscous  forces 
in  the  atmosphere  produce  heat  and  momentum  gradients.  These  gradients  cause 
nonequilibrium  conditions  where  energy  is  transported  and  dissipated  through 
turbulent  motions.  Once  created,  turbulent  eddies  gradually  break  down  into 
smaller  eddies  until  they  are  so  small  that  viscous  dissipation  into  heat 
occurs.  This  process  of  turbulence  creation  and  destruction  continues  until 
the  driving  energy  is  removed  from  the  system. 

Adequate  description  of  a  passive  substance  diffusing  into  a  turbulent 
atmosphere  requires  a  wide  variety  of  meteorological  measurements  and  algo¬ 
rithms.  The  Navier-Stokes  equations  used  to  describe  the  turbulent  atmosphere 
are  not  amenable  to  direct  solutions.  Certain  simplifying  assumptions  allow 
limited  solutions  that  can  yield  diffusion  equations  used  to  model  the  turbu¬ 
lent  atmosphere.  The  sophistication  and  accuracy  of  the  models  vary  with  the 
amount  of  detailed  turbulence  measurements  or  calculations  available  and  the 
degree  of  realism  of  the  algorithms  used. 

Several  forms  of  information  can  be  used  in  turbulence  models.  The  most 
direct  forn.  is  turbulence  measurements.  However,  these  may  not  be  available 
because  of  time,  cost,  or  logistic  problems.  Also,  most  models  lack  the 
sophistication  required  to  handle  detailed  turbulence  data.  Turbulence  can 
also  be  described  as  a  function  of  stability,  the  best  estimates  of  which  are 
ratios  of  heat  and  momentum  fluxes  derived  from  the  principal  components  of 
the  turbulent  energy  equation.  Profile  measurements  of  wind  and  temperature 
provide  mean  wind  and  temperature  model  input  data  and  are  used  for  calcula¬ 
tion  of  gradients.  If  fluxes  of  heat  and  momentum  are  assumed  to  remain  es¬ 
sentially  constant  within  the  surface  boundary  layer,  flux/gradient  relation¬ 
ships  provide  a  direct  means  of  calculating  stability  from  profile  measure¬ 
ments.  The  Obukhov  Length  is  an  excellent  stability  parameter  for  use  in 
turbulence  models.  It  is  derived  from  the  turbulent  kinetic  energy  equation 
with  a  minimum  of  qualifying  assumptions  and  can  be  reliably  obtained  from 
high  quality  wind  and  temperature  profile  measurements. 

An  adequate  description  of  atmospheric  turbulence  will  require  a  variety 
of  meteorological  instruments.  Sets  of  pyranometers  and  pyrgeometers  are 
needed  for  measurement  of  fluxes  of  long  and  shortwave  radiation  through  a 
surface.  Scintillometer  measurements  of  refractive  index  structure  near  the 
ground  will  supply  information  on  eddy  activity  in  the  viscous  dissipation 
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range.  Micrometeorological  anemometer  cups  and  vanes  and  precision  thermo¬ 
meters,  at  Togarithmically-spaced  height  Intervals  from  2  to  16  meters,  are 
the  most  reasonable  source  of  profiles  of  mean  wind  and  teiperature.  Instru¬ 
ments  capabl ?  of  measuring  high-frequency  fluctuations  could  provide  data  for 
turbulence  spectrum  analysis.  Correct  combinations  of  these  instruments  could 
provide  data  on  the  fluxes  of  heat,  momentum,  and  moisture.  Acoustic  sounders 
can  provide  in-depth  measurements  of  wind,  estimates  of  bulk  turbulence,  and 
inforiiiation  )n  the  depth  of  convection  or  mixing  layers.  Scintillometers  can 
provide  the  )0ssibility  of  calculating  path-averaged  wind  and  divergence.  An 
expanded,  improved  mesonet  grid  and  thorough  modernization  of  the  weather  sta¬ 
tion  would  provide  the  meso-  and  synoptic-scale  information  needed  to  describe 
the  setting  in  which  turbulence  events  occur. 

^  The  information  will  be  most  useful  if  the  instruments  are  deployed  and 
the  "data  collected,  analyzed,  and  entered  into  the  model  by  people  with  a  firm 
understanding  of  the  assumptions  used  and  knowledge  of  the  limitations  of  the 
site  and  conditions  of  data  collection.  When  the  data  are  collected  at  sites 
approximately  homogeneous  in  space  and  with  meteorological  conditions  station¬ 
ary  in  time,  turbulence  can  be  described  as  a  function  of  stability.  Signif¬ 
icant  departures  from  these  conditions  vastly  complicate  data  interpretation. 
Stability  parameters  can  be  used  as  test  go/no-go  criteria,  for  test  data 
quality  con.rol,  and  for  modeling.  If  data  from  the  instruments  discussed 
above  were  available,  in-depth  posttest  analysis  would  be  possible. 


1.  INTRODUCTION; 


Diffusion  meteorologists  agree  that  behavior  of  a  plume  released  into  the 
atmosphere  is  controlled  by  meteorological  influences.  Unless  the  release  is 
accompanied  by  the  emission  of  a  tremendous  amount  of  heat,  initial  source 
characteristics  are  quickly  overwhelmed  and  plume  behavior  becomes  a  function 
of  the  motions  in  the  atmosphere  into  which  the  plume  is  released.  The  mean 
wind  field  transports  the  released  material  dcwnrange  and  turbulent  motions 
cause  horizontal  and  vertical  spreading  of  the  material.  Transport  is  a  func¬ 
tion  of  mean  wind  fields  that  are  readily  measured.  Diffusion  is  a  function 
of  turbulent  eddies  that  a'"e  much  more  subtle  and  difficult  to  quantify. 

A  basic  equation  describing  concentraticn  (  x,  mg/tn^)  of  an  inert,  non¬ 
depositing  material,  diffusing  at  distance  x  (meters)  downwind  from  a  continu¬ 
ous  ground  level  gaseous  point  source,  over  a  surface  homogeneous  in  space, 
into  an  atmosphere  stationary  in  ti.ne  is 

y  ^  Q/i^  ,  (1.1) 


where  Q  is  the  source  strength  in  mg/sec ^  Oy  and  07  are  the  crosswind  and 
vertical  plume  standard  deviations  and  u  is  the  transport  wind  speed.  The 
in  Equation  1.1  is  proportional  to  the  volume  occupied  by  the  dif¬ 
fusing  material,  assuming  normal  or  Gaussian  distribution.  Derivation  of  this 
result  is  performed  in  Section  4  of  this  document.  The  double  overbar  on 
transport  wind  speed  indicates  that  it  is  a  mean  wind  speed  averaged  over 
travel  time  and  depth  of  the  cloud.  Equation  1.1  quickly  becomes  complicated 
as  different  source  configurations  or  chemical  reactions  occur,  or  as  atmo¬ 
spheric  characteristics  change. 

Pasquill  (1976)  discusses  Oy,  the  horizontal  component  of  cloud  growth, 
as  a  function  of  downwind  travel  distance  and  measured  standard  deviation  of 
the  horizontal  wind  angle  (^a).  expressed  by  equation  1.2. 


Oy  =  9  X  f(x)  (1.2) 


Equation  1.2  applies  to  a  stationary  w'nd  field  free  of  horizontal  and 
vertical  wind  shears.  These  conditions  can  be  approximated  within  the  First 
several  tens  of  meters  above  ground  level  and  over  travel  times  of  several 
tens  of  minutes.  The  function  f(x)  assumes  a  value  of  unity  near  the  source 
and  decreases  as  a  function  of  travel  distance  and  meteorological  variables. 

The  in  Equation  1.2  is  the  horizontal  wind  angle  standard  deviation, 
a  bulk  turbulence  parameter  obtained  from  turbulence  measurements.  It  should 
represent  the  integral  of  horizontal  turbulent  fluctuations  over  the  entire 
spectrum  of  eddy  sizes  for  given  time.  In  oractice,  the  measurement  is 
limited  on  the  low  frequency  (large  eddy)  side  by  sampling  duration  and  is 
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limited  on  the  high  frequency  (small  eddy)  side  by  instrunent  response  charac¬ 
teristics  anl  data  sampling  rate.  Other  factors  affecting  the  magnitude  of 
are  adequ  icy  of  instrument  response  characteristics  and  instrument 
exposure. 

Algorithis  comparable  to  Equation  1.2  are  available  for  describing  o^, 
the  vertical  measure  of  cloud  growth,  as  a  function  of  x.  The  most  useful  of 
these  algorithms  include  the  description  of  as  a  function  of  the  vertical 
component  of  turbulence,  represented  as  the  vertical  wind  angle  standard  devi¬ 
ation,  0^. 

The  scales  of  motion  that  consist  of  eddies  somewhat  smaller  than  the 
diffusing  cloud  contribute  most  significantly  to  diffusion  through  the  oa 
and  Og  terms.  Scales  of  motion  larger  than  the  cloud  tend  to  move  the  cloud 
as  a  whole.  Turbulent  eddy  activity  near  the  earth's  surface  is  the  sum  of  a 
dynamic  or  mechanical  component  generated  by  wind  blowing  over  roughness  ele¬ 
ments  (e.g.,  trees)  and  a  convective  component  caused  by  buoyancy.  The  neu¬ 
tral  or  adiabatic  atmosphere  contains  only  the  mechanical  turbulence  component 
where  eddy  activity  depends  on  the  interaction  of  wind  with  surface  roughness. 
Diabatic  atmospheres  also  include  buoyant  low-frequency  eddies,  in  the  un¬ 
stable  case,  and  an  assortment  of  high-frequency  eddies  superimposed  on  inter¬ 
mittent  low-frequency  waves  in  the  stable  case.  Algorithms  pertaining  to  the 
neutral  through  unstable  diabatic  atmosphere  near  the  surface  will  be  primar¬ 
ily  addressed  because  these  are  the  conditions  under  which  most  current  test¬ 
ing  is  done  at  US  Army  Dugway  Proving  Ground  (DPG),  Utah. 

Just  as  diffusion  calculations  require  specification  of  the  turbulence 
field,  turbulence  can  either  be  measured  or  estimated  from  a  similarity 
theory.  Surface  layer  similarity  theory  is  a  collection  of  dimensional  and 
physical  arguments  that  permit  description  of  turbulence  as  functions  of  a 
stability  parameter.  Stability  parameters  are  combinations  of  pertinent 
variables  used  to  describe  the  degree  of  departure  from  the  adiabatic  condi¬ 
tion.  The  usefulness  of  stability  parameters  is  a  function  of  the  accuracy 
with  which  they  can  be  calculated  and  how  well  they  represent  actual 
condi tions. 

This  document  describes  some  elementary  relationships  between  atmospheric 
diffusion,  turbulence,  and  stability  parameters  derived  from  similarity  argu¬ 
ments.  The  purpose  is  to  examine  these  relationships,  determine  the  measure¬ 
ments  required  for  an  adequate  description  of  the  atmosphere  for  diffusion 
applications,  and  suggest  the  types  of  equipment  that  can  be  used  to  obtain 
these  measurements. 

2.  A  DEFINITION  OF  TURBULENCE: 

The  coordinate  system  shown  in  Figure  2.1  is  used  in  the  following  discus¬ 
sion:  The  X  axis  is  positive  to  the  right,  the  y  axis  is  positive  towards  the 
top  of  the  page,  and  the  z  direction  is  a  perpendicular  directed  out  of  the 
page.  Components  of  the  wind  (u,  v,  and  w)  are  defined  along  the  x,  y,  and  z 
axes,  respectively.  The  mean  wind  flow  is  assumed  to  be  to  the  right  along 
the  X  axis.  Vector  quantities  are  indicated  with  an  arrow  overbar  ( -^ ) . 


The  stress  tensor  or  shear  stress  (t)  is  defined  as  a  force  per  unit  area 
applied  in  normal  or  tangential  components  to  a  surface.  Figure  2.1  shows  a 
volume  with  stress  tensors  applied  to  an  xy  plane  (constant  z  surface).  For  a 
viscous  fluid,  there  are  nine  stress  components,  three  on  each  plane.  Tensors 
T2y  and  Tjjx  operate  in  the  xy  plane,  while  T22  is  normal  (90°)  to  the  xy 
plane.  The  first  stress  tensor  subscript  indicates  the  face  on  which  the 
stress  acts  and  the  second  subscript  indicates  the  direction  of  stress  action. 
With  inviscid  or  frictionless  fluids,  the  tangential  stress  components  are 
zero.  Then  the  normal  stress  tensor,  acting  in  tension,  is  opposed  by  pres¬ 
sure  forces  only.  By  convention,  tension  is  designated  as  positive  and  pres¬ 
sure  (P22)  as  negative.  For  a  Newtonian  viscous  fluid,  which  the  atmosphere 
approximates,  momentum  transfer  caused  by  an  applied  tangential  stress  is 
proportional  to  the  rate  of  strain  (or  momentum  gradient).  For  a  Newtonian 
viscous  fluid  moving  horizontally  above  a  stationary  boundary,  the  stress 
^■'^zxo*  ■^zyo)  is  proportional  to  the  rate  of  strain  (3u/3z  or  3v/3z)  es 
expressed  in  Equation  2.1. 


Figure  2.1.  Coordinate  System. 


Tzxo  =  ^  3j/37.,  T^yo  =  u  3v/3z  (2.1) 


where  u  is  a  constant  of  proportionality  known  as  the  coefficient  of  viscos¬ 
ity.  The  subscript  o  on  r and  ^zyo  denotes  the  fact  that  these  tangen¬ 
tial  viscous  terms  arise  because  of  the  mctions  of  a  viscous  fluid  above  a 
stationary  surface.  In  the  atmosphere,  viscous  terms  are  small  except  near  a 
surface  where  the  velocity  gradient  is  large  or  at  the  high  frequency  end  of 
the  turbulence  spectrum.  The  x  component  of  viscous  force  per  unit  mass  on  a 
volume  of  air  is  then  represented  as  in  Equation  2.2. 


Fx  =  (1/p)  3t2^q/hz  =  (u/p)  3^u/3z^. 


(2.2) 
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The  stress  teems  described  above  are  used  in  Navier-Stokes  equations  which 
express  Newton's  Second  Law  for  motion  in  a  Newtonian  viscous  fluid.  For  the 
u  component  of  an  incompressible  viscous  fluid,  the  instantaneous  equation  of 
motion  is 

du/dt  =  (-l/p)8P/9x  +  fc  +  (u/p)[3^u/ax^  +  3^u/3y^  +  3^u/3z^].  (2.3) 


Equation  2.3  presents  acceleration  of  the  u  component  of  the  wind  as  the 
resultant  of  pressure  gradient  force,  Coriolis  force,  and  viscous  forces.  All 
forces  in  Equation  2.3  are  per  unit  mass.  The  full  vector  form  of  the  equa¬ 
tion  of  motion  is  used  in  Section  5  for  deriving  the  turbulent  kinetic  energy 
equation.  Normally,  the  rate  of  change  in  vertical  wind  speed  gradients 
greatly  exceeds  horizontal  accelerations,  that  is 

3^u/3z^  >>  3^u/3x^,  3^u/3y^  (2.4) 


and  these  smaller  terms  in  Equation  2.3  can  be  neglected.  Similar  equations 
can  be  written  for  the  v  and  w  wind  components. 

In  addition  to  the  equation  of  motion,  a  statement  of  continuity  is 
required  to  describe  the  action  of  turbulence.  The  equation  of  continuity  is 
most  easily  visualized  by  examination  of  flow  tnrough  a  unit  volume  (Figure 
2.2). 


Z 


[pu  +  (3(pu)/3x) 


4x]  AyAz 


Figur,.  "  Flow  Through  a  Unit  Volume. 


In  Figure  1.7  :onsider  mass  inflow  (pu)  multiplied  by  unit  area  AyAz  and  mass 
outflow  as  p  +  ('(ou)/3x)ax  where  (3(pu)/3x)Ax  expresses  changes  that  take 
place  while  'e  fluid  passes  through  the  volume.  Then,  the  net  rate  of  nress 
inflow  per  un  t  volutne  equals  the  local  rate  of  change  of  density,  as  given  in 
the  three-dimensional  mass  divergence  equation  (continuity  equation). 


-3p/3t  =  9(pu)/9x  +  3(pv)/3y  +  3(pw)/3z  •=  7*pV. 


(2.5) 
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Equation  2.13,  neglecting  the  total  di'rivative  and  viscous  terms,  then  becomes 


(l/p)3P/3x  -  f'c  =  (1/ P)[  3(  pKfn  3u/ 3x)/ 3x  + 
aCpKin  3u/ 3y )/ ^  -t  3(11(^30/ 3z)/ 32] 


(4.9) 


and  solutions  are  possible  with  the  assumption  that  is  some  known  func¬ 
tion  along  each  axis.  Similar  expressions  for  heat  can  be  written  with  the 
K-theory  relttionship 


w'  6'  =  Kh  30/  a-. 


(4.10) 


where  3  represents  potential  temperature. 


In  the  early  years  of  diffusion  work,  the  magnitude  o'‘‘  K  for  each  applica¬ 
tion  was  generally  unknown.  Because  of  this,  the  "Reynolds  Analogy"  was  fre- 
tiuently  employed.  This  analogy  assumed  equiva''ence  of  the  di f fusi vi ties  for 
momentum,  heat,  and  suspended  material,  i.e  , 


Kp  =  Kh  =  k  material 


(4.11) 


Scvera"!  probloms  with  this  analooy  are: 

a,  ‘Momentum  transfer  occurs  as  the  result  of  pressure  forces  which  do  not 
o;-?ratt  equivalently  on  fields  of  heat  or  suspended  material. 

[, .  A  material  released  into  the  air  at  a  temp.erature  different  than  the 
air  temiora tur-^  vvill  undergo  displacements  relati''0  to  ambient  eddy  motions. 

c.  Any  diffusing  material  motions  caused  by  heating,  evaporation,  radia¬ 
tion  absorption,  or  differential  motions  caused  by  inertial  or  aerodynamic  ef- 

ects  ('.ettlii  g  of  droplets,  sailing  of  spores,  etc.)  will  affect  diffusivity. 

d.  Diffusivity  gradients  may  also  be  caused  by  terrain  or  roughness  ele¬ 
ments  ,.'hich  I  ausG  shears,  vortex  siedding,  or  wake  flows. 

It  is,  therei  jre,  generally  accepted  tha  *■ 


K|.  a  Kf,  ^  K  ,-ia ferial 


(4.12) 


K  is  called  a  coeft-icient  of  eddy  diffusivity.  The  concept  of  K  arises 
from  an  assumption  that  the  mixing  of  eddies  in  the  atmosphere  is  analogous  to 
molecular  mixing.  Mixing  on  the  molecular  scale  occurs  is  molecules  travel 
over  a  characteristic  molecular  length  scale  defined  as  c,  mean  path  distance 
( t)  before  colliding  with  other  molecules.  Bradshaw  (1971)  points  out  several 
problems  with  the  analogy; 

a.  Turbulent  eddies  are  continuous  and  contiguous,  whereas  gas  molecules 
are  discrete  and  collide  only  at  intervals. 

b.  Although  molecular  free  raths  are  small  compared  to  the  dimension  of 
the  flow,  turbulent  eddies  are  hol. 

In  spite  of  these  deficiencies,  the  mixing  length  analogy  and  the  resultant 
e<i‘y  diffusivity  or  K-theory  relationships  remain  popular  because  they  provide 
the  simplest  solutions  to  diffusion  problems.  Alternative  numerical  methods 
are  beyond  the  means  and  requirements  of  most  operational  diffusion  programs. 

unfortunately,  all  the  complexities  of  turbulent  interactions  influence 
the  maqnit'ide  of  K.  Adequate  specification  of  K  is  difficult  because  it  is  a 
complex  function  of  height  above  the  surface  and  stability.  It  also  varies 
with  the  Characteristics  of  the  diffusing  substance  to  be  treated.  The  most 
familiar  forms  of  K  are  those  for  momentum  and  heat.  This  is  because  wind  and 
temperature  are  the  most  frequently  measured  atmospheric  variables.  The 
turbulent  transfer  coefficient  for  momentum  (K^)  is  defined  as  the  product 
of  a  characteristic  eddy  velocity  and  a  characteristic  eddy  length  scale,  e.g. 

Km=(w'w')Ht  (4.7) 


Similar  turbulent  transfer  coefficients  or  eddy  diffusivities  can  readily  be 
defined  for  sensible  heat  (K^,),  moisture  (K]),  or  fractional  atmospheric 
components  (K  material). 

An  example  of  K-theory  utility  is  the  "first  order  closure"  or  eddy 
diffusivity  (K)  closure  solutions  to  the  Navier-Stokes  equations  (Equation 
2.13),  where  eddy  stress  is  expressed  directly  in  terms  of  gradients ^obtained 
from  the  mean  motion  field.  Neglecting  the  viscous  terms  (  u/ p9  u/ 9z ^)  and 
assuming  that  the  total  change  with  respect  to  time  (du/dt)  is  zero.  Equation 
2.13  is  reduced  to  a  balance  between  the  pressure  gradient  forces,  the 
Coriolis  term,  and  Reynolds  stresses.  K-theory  relates  Reynolds  stresses  to 
mean  gradients  using  an  analogy  to  viscous  stresses  in  a  Newtonian  fluid. 
Equation  4.8,  for  example,  relates  u'w'  to  the  vertical  gradient  of  wind. 


The  existtnce  of  sou-'ces  or  sinks  which  cause  variability  in  Q  has  been 
difficult  to  cuantify.  Variations  in  Q  may  be  due  to  deposition,  physical,  or 
chemical  chanr.es.  A  diffusion  model  must  account  for  these  changes  to  satisfy 
continuity,  and  a  number  of  schemes  are  available  to  do  this.  For  some  trials 
a  "budget"  is  obtained  by  passing  a  diffusing  cloud  through  a  vertical  sampl¬ 
ing  array,  which  collects  a  series  of  point  samples  to  be  used  to  back  calcu¬ 
late  Q.  Such  sampling  procedures  tend  to  be  elaborate  and  involve  consider¬ 
able  expense.  Also,  cloud  growth  rates  constrain  this  procedure  to  use  near 
the  source,  /another  approach  is  to  take  advantage  of  modern  tracer  technol¬ 
ogy.  Modern  tracers  are  inert,  nontoxic  gases  that  have  low  background  con¬ 
centrations  arid  are  readily  measured  to  low  concentrations.  When  released  in 
a  known  ratio  with  other  substances,  these  tracers  can  be  used  to  evaluate  the 
presence  of  sources  or  sinks  for  the  substance  being  tested.  Remote  sensing 
techniques  such  as  lidar  scanning  may  also  help  resolve  source  strength  mea¬ 
surement  problems. 


Assuming  nondivergence  of  the  mean  component,  as  in  Equation  4.3,  the 
total  derivative  of  xwith  the  nonzero  mean  components  fro'';  Equation  4.2 
becomes 


-d  x/dt  =  -[  3x/  ^  +  11  9x/  +  V  3y/  ay  +  W  ay/  ^z] 

=  3{"u''x')/9x  +  (V'T)/^  +  atw‘l(’)/3z 


(4.4) 


At  this  point,  K-theory  or  gradient  transfer  theory  is  used  to  obtain  an 
operational  solution  for  some  models.  The  alternatives  are  cumbersome 
numerical  techniques  that  require  extensive  computer  time  and  have  not 
demonstrated  practical  utility  in  solving  diffusion  problems.  With  K-theory, 
it  is  assumed  that  the  flux  of  a  substance  across  a  fixed  surface  can  be 
expressed  in  terms  of  the  mean  gradients  of  that  substance  normal  to  the 
surface,  with  a  diffusivity  term  acting  as  a  constant  of  proportionality.  For 
example,  the  first  term  on  the  right  in  Equation  4.4,  (u'x')  can  be  expressed 
by 


u'x'  -  -  Kx  m  3x,  (4.5) 


Where  Kx  is  the  diffusivity  coefficient  along  the  x-axis.  ' 

Analogous  expressions  are  readily  found  for  the  other  components.  Replacing  ' 

the  eddy  flux  terms  with  gradient  transfer  forms,  the  right  side  of  Equation 

4.4  becomes  1 

1 

dl^dt  =  3(Kx  3x/3x)/ax  I  8(Ky  3l/3y)/3y  +  3{Kz  3l/3z)/3z  (4.6)  ' 

1 

1 

1 
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accumulated  dose,  etc.  The  remainder  of  this  section  will  discuss  only  x  in 
its  relation  to  diffusion  equations,  because  sampling  techniques  are  beyond 
the  scope  of  this  document. 


Consider  a  fractional  component  or  "pollutant"  emanating  from  a  point  or 
line  source  in  the  atmosphere.  At  the  source,  the  concentration  approaches 
infinity  (xo  ”)  and  at  infinite  distance  from  the  source,  the  concentra¬ 
tion  approaches  zero  (x®  +0).  At  intermediate  distances,  the  gradients 
3x/3x,  3x/3y,  and  3x/3z  are  finite  in  the  x,  y,  and  z  directions.  Within 
these  limits,  a  statement  of  continuity  for  an  incompressible  fluid  analogous 
to  Equation  2.5  pertains. 


-3x/3t  =  3{ux)/3x  +  3(vx)/3y  +  3{wx)/3z  =  V*xV  =  V»Vx  +  xV*V.  (4.2) 


The  X  in  Equation  4.2  is  composed  of  a  mean  (x)  component  representative  of 
the  average  concentration  in  the  volume  of  interest,  and  a  deviation  (x')  due 
to  eddies  carrying  different  concentrations  as  they  move  through  a  measured 
volume.  The  product  of  the  fluctuating  (u',  v',  w')  component  with  x'  repre¬ 
sents  a  flux  through  the  volume.  As  before,  average  values  of  fluctuating 
components  are  zero,  but  correlated  products  of  fluctuating  components  are 
nonzero. 

At  this  point  the  first  set  in  a  series  of  assumptions  are  made: 

a.  The  assumption  of  nondivergence  is  applied  to  the  mean  concentration 
component. 

x{3u/3x  +  3v/3y  +  3w/3z)  =  =  0.  (4.3) 

I 

,’w 

b.  Density  is  assumed  to  be  constant.  Variations  in  x  are  then  essen¬ 
tially  caused  by  varying  amounts  of  the  fractional  component,  not  by  changes 
in  air  density. 

c.  It  is  assumed  that  there  are  no  additional  sources  or  sinks  for  the 
fractional  component,  i.e.,  Q  remains  constant. 

Measurements  should  be  made  in  every  diffusion  test  to  check  the  degree  of 
conformity  of  actual  conditions  with  these  assumptions.  Divergence  computa¬ 
tions  from  path-averaged  wind  measurements  could  verify  the  first  assumption. 

If  possible,  the  divergence  computation  should  be  made  on  length  and  time 
scales  similar  to  the  scale  of  the  diffusion  test.  Although  the  magnitude  of 
divergence  does  not  directly  enter  the  diffusion  equations  used  to  model  trial 
results,  the  presence  of  significant  divergence  could  explain  some  of  the 
deviations  of  the  actual  measurements  from  expected  results. 

Density  is  easily  monitored  as  a  function  of  temperature  and  pressure.. 

Density  variations  at  a  specific  site  during  a  test  series  are  usually  incon¬ 
sequential.  However,  density  differences  must  be  considered  when  test  results 
from  different  sites  or  under  significantly  different  meteorological  condi¬ 
tions  are  compared. 
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Uj/Ui  =  (Z2/Z1)'" 


(3.14) 


where  m  is  soiie  fractional  exponent.  Exponent  m  assumes  a  value  near  0.14  for 
the  neutral  case,  but  varies  with  stability,  site  roughness,  and  height  of 
measurement  above  the  surface.  Exponent  m  is  difficult  to  specify  accurately 
because  its  magnitude  is  influenced  by  many  different  factors.  This  is 
unfortunate  because  small  variations  in  m  cause  large  changes  in  U2/U1.  In 
spite  of  these  limitations,  power  law  formulations  are  popular  because  they 
provide  the  simplest  solutions  to  diffusion  equation  wind  profiles.  An 
accurate  specification  of  m  for  diffusion  tests  requires  continuous,  repre¬ 
sentative  readings  of  wind  profiles  through  the  depth  of  a  diffusing  cloud. 
Remote  wind  sounding  techniques  are  wall  suited  to  provide  these  data  within 
the  lower  troposphere. 

4.  K-THEORY  AND  GAUSSIAN  DIFFUSION: 

With  a  basic  definition  of  turbulence  and  two  approaches  to  SBL  wind  pro¬ 
files  at  hand,  it  is  now  appropriate  to  examine  some  basic  diffusion  relation¬ 
ships.  This  is  not  a  critical  examination  of  diffusion  models,  for  which  many 
references  exist.  This  is  a  discussion  of  some  of  the  assumptions  used,  but 
frequently  not  noted,  in  the  process  of  developing  models  or  designing  diffu¬ 
sion  test  programs.  Suggestions  are  also  made  concerning  measurements  neces¬ 
sary  to  verify  assumptions  and  conduct  satisfactory  diffusion  tests. 

The  quantity  of  airborne  material  is  frequently  described  or  measured  as 
a  fractional  component  of  the  air  in  units  of  parts  per  million  or  similar 

units  easily  converted  to  a  mixing  ratio  (R).  For  example,  the  National 

Institute  of  Occupational  Safety  and  Health/Occupational  Safely  and  Health  Act 

(NIOSH/OSHA)  Pocket  Guide  to  Chemical  Hazards  (Mackison  and  Stricoff,  1980) 
describes  permissible  exposure  limits  for  pollutants  in  this  form.  However, 
diffusion  equations  require  introduction  of  fractional  components  in  terms  of 
a  concentration  (x)  with  units  of  (milli)grams  per  cubic  meter,  or  related 
units.  The  conversion  from  R  to  x  is  achieved  by  multiplying  the  number  of 
(milli)grams  of  fractional  component  per  gram  of  air  by  air  density  [in 
(millilgrams  per  cubic  meter]. 


X  =  Rp  • 


(4.1) 


R  is  usually  measured  by  instruments  which  must  sample  over  finite  time 
intervals.  When  sampling  duration  must  be  considered,  the  units  of  measure 
become  a  time-averaged  concentration  or  a  dosage  (dd)  with  units  of  (milli)- 
gram*minutes  per  cubic  meter.  Also,  a  measure  of  the  total  accumulation  of  a 
fractional  component  over  a  period  of  time  per  unit  area  is  often  desirable. 
The  lungs,  for  example,  act  as  collecting  devices,  accumulating  inhaled  compo¬ 
nents  over  time.  Many  deposition  or  impaction  type  samplers  operate  in  this 
fashion,  but  with  varying  efficiencies,  yielding  a  measure  of  total  accumu¬ 
lated  dose  (D)  with  units  of  (milli  )grams.  Total  dose  samplers  are  usually 
fairly  simple  collection  devices,  but  their  efficiencies  are  complicated  by 
factors  such  as  wind  speed,  angle  of  exposure  to  the  wind,  particle  size. 
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The  diabati ;  Influence  function  accounts  for  nonlinearities  in  the  logarithmic 
wind  profile,  caused  by  thermal  stratification.  For  neutral  conditions  4ini  = 
1.0  and  ij)  =  0.0.  Table  3.1  gives  sample  values  of  (t)^  and  as  functions  of 
the  stability  parameter  z/L  (defined  in  Section  6). 


Table  3.1.  Dimensionless  Shear 
{((im)  ahd  the  Diabatic  Influ¬ 
ence  Function  (i|/)  as  Functions 


of  the 

Stability 

Parameter  z/L 

for  the 

Unstable 

Case. 

z/L 

I’m 

0.0 

-0.1 

0.795 

-0.2 

0.707 

0.442 

-0.5 

0.586 

-1.0 

0.500 

1.084 

-5.0 

0.339 

2.023 

The  bouhdaries  of  the  SBL,  where  stress  remains  nearly  constant  and  the 
logarithmic  law  applies,  can  be  defined.  The  lower  boundary  is  Zq,  the  top 
of  the  viscous  sublayer  where  viscous  shear  becomes  dominant  at  the  expense  of 
turbulent  shear  stress.  Garratt  (1980)  finds  that  viscous  effects  distort  the 
wind  profile  for  heights  in  the  range  l<z/zo<100  and  suggests  a  method  to 
compensate  for  this  effect.  The  upper  boundary  is  the  top  of  the  SBL.  The 
top  of  the  SBL  is  somewhat  arbitrarily  defined  as  the  height  at  which  t 
departs  from  tq  by  more  than  20  percent.  For  a  representative  u*  of  50  cm 
sec"^  and  density  held  constant  at  0.0012g  cm“^ ,  tq  has  a  value  of 
3.0  g  cm"^  sec"^  .  Using  Equation  2.13  and  ignoring  the  viscous  term,  non¬ 
accelerating  components  are  in  balance  when  (-l/p)9P/9x  and  f^  are  equal  and 
of  opposite  sign.  Then,  accelerations  (du/dt)  are  due  to  spatial  changes  in 
Reynolds  stress  terms.  If  changes  in  the  horizontal  stress  components  are 
neglected  (as  assumed  earlier),  accelerations  are  balanced  by  changes  in  stress 
with  height.  A  typical  wind  speed  change  of  1.0  m  sec"^  in  10  minutes 
produces  a  du/dt  of  0.167  cm  sec~^ .  At  this  rate,  stress  remains  within  20 
percent  of  tq  for  30  meters.  It  can  be  generally  assumed  that  constant  stress 
and  the  logarithmic  wind  profile  are  useful  at  DPG  in  the  layer  ranging  from  a 
few  centimeters  to  a  few  tens  of  meters  above  the  ground. In  some  cases  the 
logarithmic  wind  profile  ray  be  representative  up  to  100  meters  above  ground 
level,  well  beyond  the  top  of  the  SBL.  The  logarithmic  wind  profile  and 
attendant  similarity  assumptions  can  therefore  apply  in  the  layer  of  interest 
for  the  majority  of  diffusion  testing  at  DPG. 

The  logarithmic  wind  profile  (Equation  3.12)  should  be  contrasted  with  the 
empirical  power  law  wind  profile. 


a 


E 
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Because  ot  the  influence  of  roughness  length  on  the  wind  profile,  test 
sites  for  turoulence  moasurements  must  be  carefully  chosen.  Variations  in 
roughness  witn  wind  fetch  greatly  coiriplicate  interpretation  of  turbulence 
measurements.  Also,  introduction  of  nan-made  obstacles  should  be  kept  to  a 
minimum  at  turbulence  measurement  sites.  For  sites  with  trees,  large  bushes, 
or  similar  impediments  to  flow,  a  correction  factor,  displacement  height  (d), 
must  be  introduced.  Displacement  height  is  a  function  of  the  height  and  the 
spacing  of  the  flow  impediments.  Like  Zq,  d  is  a  characteristic  of  flow  at 
a  specific  site,  and  may  suffer  from  aerodynamic  effects  similar  to  Zq  in 
strong  winds.  The  least  complicated  test  sites  for  turbulence  measurements 
contain  no  obstacles  to  cause  a  requirement  for  d.  At  DPG,  bushes  are  widely 
spaced  and  generally  less  than  half  a  meter  tall.  Therefore,  d  is  considered 
inconsequential  because  of  the  general  absence  of  large  or  densely  packed 
obstacles  to  the  flow.  The  full  form  of  the  logarithmic  wind  profile  equation 
for  an  adiabatic  atmosphere  is  given  by  Equation  3.10,  with  u(z)  represent¬ 
ing  a  mean  wind  at  height  z  (valid  for  z>(d+Zo)). 


u(2)  =  {u*/k)  ln(  (z-d-Zo)/Zf..) 


(3.10) 


•  At  well  chosen  test  sites,  Zq  and  d  are  small  enough  to  be  neglected  in  the 

numerator  ot  Equation  3.10. 

. To  this  point,  the  effects  of  stability  on  the  wind  profile  have  not  been 
discussed.  For  near-neutral  conditions,  a  plot  of  wind  speed  versus  height  is 
nearly  a  straight  line  on  a  log-linear  plot,  as  in  Figure  3.1.  However,  the 
B  wind  profiles  show  significant  departures  from  linearity  for  diabatic  condi¬ 

tions.  During  unstable  conditions,  the  wind  profile  is  expected  to  increase 
more  slowly  than  in  the  neutral  case  and,  in  stable  conditions,  the  wind  pro- 
file  increases  more  rapidly  than  in  the  neutral  case  (Figure  3.1).  To  account 
for  departures  from  linearity,  a  new  variable,  <)^,  the  dimensionless  wind 
^  shear,  is  introduced  to  Equation  3.9  so  that 


3u/3z  =  u*i^/kz,  (3.11) 


:• 


and  a  diabatic  influence  function  (  i|j)  derived  as  a  function  of  by  Paulson 
(1970)  is  introduced  to  the  logarithmic  wind  profile,  as  shown  in  Equation 
3.12. 


u[z)  -  (u*-/k)  ln((z/zo)-'i') 


(3.12) 


Equation  3.13  presents  the  functional  relationship  of  to 

^=2  ln((l  +  <l|^^)/2)  +ln((l+  ^^)/2)  -2tan-^(|t^W  tt/2 


(3.13) 


(u2-ui)/u*  =  1n(z2/zi  )/k. 


(3.8) 


where  k  is  an  empirical  constant  known  as  von  Karman's  constant.  Wind  profile 
studies  showing  a  linear  relationship  between  u  and  In  z  under  adiabatic 
conditions  tend  to  confirm  the  validity  of  Equation  3.8.  The  derivative  of 
Equation  3.8  over  an  infinitesimal  height  yields  the  basis  of  a  fundamental 
flux-gradient  relationship,  given  in  Equation  3.9. 


3u/3z  =  u*/kz  (3.9) 


If  Zi  in  Equation  3.8  is  chosen  sufficiently  close  to  the  ground,  u^  is  at 
the  top  of  the  viscous  sublayer  and  is  free  of  Reynolds  stress  momentum  trans¬ 
fers  from  above.  This  sublayer  height  is  determined  by  the  dynamic  roughness 
of  the  underlying  surface.  The  dynamic  roughness  height  (zq),  according  to 
Monin  and  Obukhov  (1954),  defines  the  characteristic  scale  of  the  underlying 
surface's  micro  nonuniformities.  The  magnitude  of  Zq  is  determined  by  the 
properties  of  the  underlying  surface  alone  and  should  be  independent  of  wind 
speed.  In  practice,  Zq  changes  somewhat  as  a  function  of  wind  speed  because 
as  wind  speed  increases,  plants  or  other  roughness  elements  may  assume  more 
aerodynamically  efficient  shapes  to  reduce  drag. 

Determination  of  Zq  for  a  given  site  is  done  experimentally,  with  care¬ 
ful  measurements  near  the  surface  by  using  fast  response  instruments,  such  as 
hot  wire  anemometers.  For  very  flat  surfaces,  Zq  can  be  less  than  one  cen¬ 
timeter.  Conversely,  Zq  can  approach  1  meter  in  built-up  urban  areas.  For 
near  neutral  wind  profiles,  Zg  may  be  obtained  at  the  intersection  of  the 
linear  wind  profile  at  u  =  0  on  a  plot  of  u  versus  In  z  (Figure  3.1).  Esti¬ 
mates  of  Zq  are  obtained  from  the  averages  of  repeated  measurements  using 
this  procedure.  Rough  estimates  of  Zq  can  also  be  obtained  from  the  general 
characteristics  of  the  site,  as  demonstrated  by  Hogstrom  and  Hogstrom  (1978). 
Analysis  of  wind  fields  at  DPG  yield  Zg  values  ranging  from  2  to  4  cm. 


100.0 
ln(z)  10.0 
1.0 
0.1 


Figure  3.1.  Wind  Speed  Versus  Height  (logarithmic  scale). 


The  SBL,  extending  from  a  few  centimeters  to  several  tens  of  meters  above  the 
surface,  is  of  great  interest  in  turbulence  and  diffusion  studies  because  this 
is  the  layer  in  which  most  hjman  activities  occur.  Within  the  SBL  the  ratio 
t/p  is  nearly  constant.  This  ratio  has  units  of  m^sec"  ,  the  square  root  of 
which  yields  units  of  velocity.  Consequently,  the  square  root  of  t/p  is  de¬ 
fined  as  the  friction  velocity  (u*).  In  the  SBL, 

=  t/  p  H  u*.  (3.3) 


Friction  velocity  is  the  characteristic  velocity  within  the  SBL.  It  accounts 
for  the  effects  of  large  scale  pressure  gradients  and  surface  roughness 
effects  on  SBL  winds  (see  Tennekes,  1982).  Friction  velocity  is  used  to  form 
dimensionless  ratios  with  other  SBL  winds.  The  nondimensional  SBL  winds  can 
then  be  conveniently  described  as  a  function  of  height,  surface  roughness,  and 
stability. 

Monin  and  Obukhov  (1954)  used  similarity  between  laboratory  flow  and  flow 
within  the  SBL  to  obtain  a  logarithmic  form  of  the  wind  profile.  Within  the 
SBL,  the  mean  wind  speed  difference  at  two  heights  is  related  to  the  ratio  of 
the  heights  as  expressed  in  Equation  3.4, 

(u2-ui)/u*  =  f  (za/zi)  (3.4) 


where  f  represents  a  function.  To  find  the  nature  of  f,  three  wind  speed  and 
height  combinations  were  chosen  such  that  Z3>2  2>zi  and  U3>U2>oi. 

Because 


and 


U3  -  01=  (U3  -  U2)  +  (U2  -  ui) 


(3.5) 


z  3/z  1  =  (z  3/Z2)  ’(zz/z  1), 


(3.6) 


the  functional  form  that  satisfies  Equations  3.5  and  3.6  is  the  logarithmic 
equality  presented  in  generalized  form  using  Equation  3.7. 


ln(  C3/C2)  •ln(  52/ 5i)  =  ln(  C3/?i)  =  lnC3  *  InEi  = 
(InCs-  InCz)  +  (InCz  -  InCi) 


(3.7) 


Therefore,  the  required  functional  relationship  between  the  nondimensional ized 
wind  of  Equation  3.4  and  the  height  ratio  is. 
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L.  F.  Richardson  is  said  to  have  stated  it  as  follows; 


Great  whorls  have  little  whorls 
That  feed  on  their  velocity... 

And  the  little  whorls  have  lesser  whorls. 
And  so  on  to  viscosity. 


3.  THE  LOGARITHMIC  AND  POWER  LAW  WIND  PROFILES: 

Two  empirical  equation  forms,  the  logarithmic  and  power  law  profiles, 
describe  the  change  in  wind  speed  with  height.  The  logarithmic  profile  origi¬ 
nated  with  laboratory  fluid  flow  measurements  and  has  become  part  of  the  sys¬ 
tem  of  equations  used  for  similarity  theory  descriptions  of  the  turbulent 
atmosphere.  Power  law  profiles  came  from  attempts  to  fit  vertical  wind  speed 
measurements  with  a  simple  exponential  equation.  The  simpler  power  law  forms 
are  more  mathematically  tractable  for  use  with  diffusion  equations  than  are 
logarithmic  profiles. 

To  overcome  the  impasse  presented  by  intractable  Navier-Stokes  equations, 
fluid  dynamicists  used  laboratory  experiments  for  engineering  solutions  to 
fluid  dynamics  problems.  A  number  of  very  useful  semi -empirical  relationships 
were  developed  from  these  experiments.  One  is  the  Reynolds  Number  (RE),  a 
dimensionless  ratio  of  accelerations  to  viscous  stress  gradients.  When  RE  is 
large,  flow  is  considered  to  be  fully  turbulent  regardless  of  the  scale  of  the 
experiment.  With  a  large  RE,  the  principle  of  dynamic  similitude  can  be  used 
to  relate  laboratory  findings  to  turbulent  processes  in  the  atmosphere. 

Dynamic  similitude  considerations  led  to  the  application  of  logarithmic  fluid 
velocity  profiles  observed  in  the  laboratory  to  wind  flow  in  the  friction 
layer  near  the  ground. 

Several  unique  flow  layers  have  been  found  in  laboratory  flow  simulations. 
In  pipe  flow,  there  is  a  thin  layer  (the  viscous  sublayer)  close  to  the  wall 
where  the  influence  of  viscosity  is  great  and  the  viscous  shear  stress 
(u  3u/3z)  is  much  larger  than  the  Reynolds  shear  stresses.  The  viscous 
sublayer  is  often  called  the  laminar  sublayer  due  to  absence  of  Reynolds 
stress.  Stress  at  the  wall  (tq)  is  given  by  Equation  3.1,  where  u  is  now 
defined  as  a  mean  wind  along  the  direction  of  flow. 


Tq  =  P  3u/3z 


(3.1) 


Outside  the  viscous  sublayer  in  pipe  flow,  stress  remains  constant  at  its 
wall  value  (t^tq),  but  Reynolds  stresses  are  predominant  and  viscous  effects 
are  small.  Soviet  scientists  Monin  and  Obukhov  (1954)  noted  the  similarity  of 
this  pipe  flow  phenomenon  to  observed  flows  in  the  atmosphere.  They  defined 
the  atmospheric  constant  stress  layer,  also  known  as  the  surface  boundary 
layer  (SBL),  as  a  layer  several  tens  of  meters  in  depth  in  which 


t=tq  =  -pu'w'  =  constant. 


(3.2) 
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When  winds  are  divided  into  mean  (u)  and  fluctuating  (i  ')  components,  the 
cssumption  is  made  that  the  sampling  time  is  sufficiently  hort  so  variations 
^n  the  mean  motion  are  negligible,  and  sufficiently  long  sfj  that  the  average 
value  of  u'  approaches  zero.  In  other  words,  a  "spectral  'jap"  is  assumed  to 
exist  between  the  mean  and  turbulent  components  of  flow,  'n  practice,  this 
condition  can  be  approached  over  uniform  terrain  in  the  atmospheric  surface 
layer  with  sampling  times  of  10  to  60  minutes,  assuming  thrit  the  wind  field  is 
not  perturbed  by  transient  mesoscale  features.  Under  the  assumption  that  a 
spectral  gap  exists,  the  flow  is  statistically  stationary,  i.e.,  turbulence 
statistics  of  the  flow  do  not  change  with  time  and  may  be  adequately 
represented  by  time  averages  obtained  from  equipment  mounted  at  fixed 
locations  past  which  the  fluid  travels. 

To  illustrate  energy  transfer  from  mean  to  turbulent  flow,  consider  an 
average  size  eddy  subjected  to  finite  positive  3w/ 3z  in  a  fluctuating  wind 
field  containing  eddies  larger  and  smaller  than  the  eddy  in  question.  If  eddy 
volume  remains  constant,  stretching  of  the  eddy  'n  the  z  direction  by  the 
notions  of  larger  eddies  will  be  accompanied  by  contraction  along  the  horizon¬ 
tal  axis.  Conservation  of  angular  ntomentum  requires  the  product  of  vorticity 
and  radius  to  remain  constant  in  the  absence  of  v.scous  dissipation.  For 
example,  conservation  of  angular  momentum  is  the  cause  of  the  increase  in  spin 
of  an  ice  skater  who  brings  his/her  arms  in  close  to  his/her  body.  Conse¬ 
quently,  the  response  of  an  eddy  to  3w/  3z  stretching  is  an  increase  in  the 
horizontal  component  of  eddy  vorticity  as  the  horizontal  dimensions  of  the 
eddy  decrease.  Vertical  stretching  therefore  causes  an  anomaly  in  the  hori¬ 
zontal  wind  field.  As  the  larger  eddies  stretch  the  sample  eddy,  smaller 
eddies  are  being  stretched  by  the  sample  eddy,  diluting  its  vorticity.  The 
length  scale  of  eddies  decreases  with  stretching  and  mixing  at  each  stage, 
establishing  a  cascade  of  energy  to  smaller  eddy  sizes. 

As  energy  progressively  cascades  to  smaller  length  scales,  there  is  an 
increasing  tendency  for  equiparti tion  of  eddy  energy  along  each  axis.  The 
condition  of  local  homogeneity  is  achieved  at  the  scale  where  the  eddy  statis¬ 
tical  properties  depend  only  on  spatial  separations  and  time  differences. 
Isotropy  occurs  as  the  eddy  properties  become  independent  of  axes  orienta¬ 
tion.  The  eddy  size  at  which  isotropy  occurs  is  limited  by  height  above  the 
ground.  The  probability  that  larger  eddy  sizes  will  exhibit  isotropy  in¬ 
creases  with  height  above  the  surface.  For  measurements  at  a  given  height, 
the  largest  eddy  size  for  which  isotropy  occurs  represents  the  beginning  of 
the  inertial  subrange.  The  spectrum  of  eddy  sizes  in  the  inertial  subrange  is 
a  function  of  the  rate  of  turbulent  energy  dissipation  per  unit  mass  (e).  At 
progressively  smaller  eddy  sizes,  the  inertial  subrange  merges  into  the  vis¬ 
cous  subrange,  where  the  eddy  size  spectrum  is  a  function  of  u/p  and  e.  At 
the  lower  end  of  the  viscous  subranoe,  eddy  energy  is  converted  to  the  kinetic 
energy  of  molecular  motion  (heat).  Thus  the  nonlinear  Reynolds  stress  terms 
of  the  Navier-Stokes  equations  extract  energy  from  the  mean  flow  and  cascade 
it  down  the  spectrum  of  eddy  sizes.  The  viscous  term  acts  as  a  sink  at  the 
lower  end  of  the  spectrum. 

Following  the  above  arguments,  Bradshaw  (IQ^l)  defines  turbulence  as  a 
three-dimensional  time-dependent  notion  in  which  vortex  stretching  causes 
velocity  fluctuations  to  spread  to  all  wavelengths  between  a  minimum  deter¬ 
mined  liy  viscous  forces  and  a  maximum  determined  by  the  boundary  conditions  of 
the  flow. 


A  more  useful  form  of  Equation  2.10  is  obtained  by  assuming  nondivergence  for 
the  fluctuating  components. 


3u'/3x+^'/^+3w'/32  =  o,  (2.11) 

Multiplying  Equation  2.11  by  H ' ,  and  recognizing  that 


3(u'17')/^  +  3(ir'7')/^+  3Cu'w')/3z  =  u'  Sr7  3x  +  u 
F  S77  3z  +  ir' ^7  ax  +  7'  Sr7  ^  +  w'  Sr7  32, 


9v  7  3y  + 

(2.12) 


the  variance  and  covariance  terms  on  the  left  side  of  Equation  2.12  can  be 
su!)stituted  into  the  fluctuating  part  of  Equation  2.10.  This  produces  the 
mean  equation  of  motion  in  an  atmosphere  containing  mean  and  fluctuating 
components  as  shown  in  Equation  2.13. 


du/dt  =  (-1/  p)  ap/  3x  +  fc  +  (  P/  p)  3^/  3z^  +  (1/  p)  3(-  pu'F  )  /  3x  + 
(1/p)  3(- (xr'7' )/ 3y  +  (1/p)  3(-piJ'w' )/ 3z 


(2.13) 


The  mean  terms  in  Equation  2.13  are  similar  to  those  in  Equation  2.3,  with  the 
major  difference  being  the  addition  of  mean  fluctuating  component  products. 

The  additional  terms  formed  from  the  variances  and  covariances  of  fluctuating 
wind  components  are  known  as  Reynolds  stresses.  In  the  atmosphere,  Reynolds 
stresses  are  larger  than  the  viscous  terms  of  the  Navier-Stokes  equations 
except  very  near  the  surface.  The  vertical  terms  containing  -pu'w'  are  of 
greatest  interest  because  the  rate  of  change  in  stress  in  the  x  and  y  direc¬ 
tions  is  usually  assumed  to  be  small  [i.e.,  3(u'u')/3x,  3(u'v')/^  approach 
zero].  Shear  stress  is  defined  in  Equation  2.14  for  the  u'w'  covariance 
term. 


Tzx  =  -pu'w'  =  -Fn, 


(2.14) 


The  u'w'  of  Equation  2.14  must  be  represented  by  ((u'w')^+  (v'w*^))^  in  the 
atmoshere.  Equation  2.14  describes  the  mean  transport  of  u'  per  unit  time 
through  a  horizontal  plane,  which  is  a  the  vertical  component  of  momentum  flux 
(F^)  due  to  Reynolds  stress.  Comparison  of  Equations  2.1  and  2.14  shows 
that  a  gradient  imposed  on  a  flow  containing  mea'  and  fluctuating  components 
will  cause  a  momentum  flux  to  develop.  Flux-gradient  relationships  are  devel¬ 
oped  further  in  Section  6. 
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In  Figure  2. 3. a,  3u/ 9x  +  ^v/ ^  is  negative,  and  3w/ 3z  must  be  equal  and  of 
opposite  sign  if  the  total  velocity  vector  remains  nondi vergent.  This 
condition,  where  the  horizontal  component  is  negative,  is  known  as 

convergence.  Figure  2.3.b  demonstrates  the  opposite  effect  for  divergence 
(Vf^*V  positive  and  3w/ 8z  negative).  In  addition  to  convergence/divergence, 
a  fluid  cross  can  display  rotation  about  the  origin,  or  vorticity  (Figure 
2.4). 


Figure  2.4.  Vorticity  (  3v/ 3X  -  3u/3y). 


The  divergent?  and  vorticity  concepts  introduced  above  are  used  with  the 
equations  of  motion  to  describe  the  action  of  turbulence.  First,  however, 
additional  stress  terms  caused  by  atmospheric  eddy  motions  must  be  introduced. 
Consider  an  instantaneous  wind  field  (u)  divided  into  if.  mean  (u)  and  fluctu¬ 
ating  (u')  components  such  that  the  time  average  of  individual  fluctuating 
components  is  zero  (u'  =  0),  but  the  products  u'u',  u'v',  etc.  can  be  nonzero. 
Then,  the  x  component  of  mass  inflow,  or  momentum  flux,  through  area  AyAz  in 
Figure  2.2  is 


u(Pu^Az)  =  pu  AyAz  =  p{u  +u 


)  ‘"Ay  Az 


=  p{u  ^ 


+  2uu'  +u'  )  £yAz  (2.9) 


The  mixed  products  (such  as  uu')  average  to  zero,  but  the  average  of  products 
of  the  fluctuating  terms  (u'w',  etc.)  are  nonzero  so  long  as  there  is  some 
degree  of  correlation  between  the  fluctuating  components.  These  nonzero 
products  must  be  considered  in  the  equations  of  mean  motion. 

To  show  the  contribution  of  fluctuating  variance  and  covariance  terms  to 
the  mean  Navier-Stokes  equations,  divide  the  total  derivative  du/dt  (from 

Equation  2.3)  into  its  local  (  3u/ 3t)  and  advective  (V  ’Vu)  components  and 
expand  it  into  mean  (u)  and  average  fluctuating  (u*)  components.  The  nonzero 
products  are 


du/dt  =  3u/  3t  +  LI  3u/  3x  +  V  3u/  ^  +  w  9u/  3z 
+  u'  3u73x  +  v'  3u73y  +  w’  3u73z 


(2.10) 


Because 


-  3p/  3t  =  V*oV  =  V  •  Vp  +  pV*V , 


(2.6) 


and  assuming  horizontal  changes  in  density  are  small  ( p  constant  in  spac.j  and 
time)  over  small  distances  in  the  surface  boundary  layer  (SBL),  it  is  reason¬ 
able  to  assume  that  the  fluid  is  homogeneous  and  incompressible.  Then, 


-3p/3t  =  0,  V*Vp  =  0,  pv»v  =  0. 


(2.7) 


Of  the  equations  represented  in  Equation  2.7,  the  divergence  equation  ( pV*V  = 
o)  is  the  one  of  interest.  Equation  2.8  is  the  divergence  equation  for  a  non- 
divergent  velocity  vector. 


V*V  =  3u/  3x  +  3v/  3y  +  3w/  3z  =  0 


(2.8) 


Eddy  motions  in  the  atmosphere  are  conveniently  described  with  the  aid  of 
a  "fluid  cross"  which  defines  the  area  occupied  by  a  slice  through  a  hypothet¬ 
ical  eddy.  A  two-dimensional  representation  of  the  hypothetical  eddy  under¬ 
going  horizontal  convergence  and  divergence  is  shown  in  Figures  2. 3. a  and  b. 
Arrows  in  the  figures  represent  velocity  components  along  each  axis,  depicting 
3u/  3x  and  3v/ 


Figure  2. 3. a.  Convergence  (  3u/ 3x  +  3v/ 9y  negative,  horizontal  area  decreases) 
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Figure  2.3.b.  Divergence  (  3u/ 3x  +  3v/ 3y  positive,  horizontal  area  increases) 


For  practical  purposes,  K  is  usually  specified  in  diffusion  equations  for  an 
inert  passive  scalar  of  neutral  buoyancy.  Other  adjustments  are  then  made 
within  models  to  account  for  departures  from  this  condition. 

With  the  foregoing  discussion  of  K  in  mind.  Equation  4.6  can  be  configured 
for  application  to  specific  types  of  diffusion  problems.  If  turbulence  is 
assumed  to  be  isotropic,  i.e.,  negligible  variation  in  diffusivity  with 
respect  to  spatial  orientation,  then 


Kx  -  Ky  =  Kz  =  K  (4.13) 


When  Equation  4.13  is  satisfied  and  K  is  constant,  simple  or  Fickian  diffusion 
occurs.  For  the  limiting  condition  where  V  equals  zero  and  Equation  4.13 
holds.  Equation  4.6  reduces  to  the  Laplacian  (v^)  of  concentration  multiplied 
by  K 

3x/8t  =  K(32x/3X^  +  3^x/3/  +  3^x/3Z^  )  =  KV^X-  (4.14) 


Equation  4.14  applies  to  instantaneous  point  source  releases  such  as  a  puff 
formed  by  the  air  burst  of  a  smoke  munition  under  calm  wind  conditions. 

For  the  more  general  condition  where  S  *  0.,  the  axis  of  plume  travel  is 
chosen  so  that  it  is  along  the  axis  of  the  mean  horizontal  wind.  With  the 
condition  that  the  mean  vector  wind  is  finite  only  along  the  x  axis,  the 
vector  wind  equation  is  reduced  to  Equation  4.15 

V  =  tu;  with  l<w  =  0,  (4.15) 


where  i,  j,  k  are  unit  vectors  in  the  x,  y,  and  z  directions. 
Then  Equation  4.6  becomes 


(3x/3t  +  u  3x/3x)  =  3(Kx  3x/3x)/9x  +  3 (Ky  3x/9y)/9y  +  9 (K^  9x/9z)/3z.  (4.16) 


If  a  further  assumption  is  made  thatsteady  state  conditions  exist,  the  first 
term  on  the  left  of  Equation  4.16  (3y/3t}  is  negligible.  For  a  continuous 
point  source  (such  as  a  plume  from  a  smoke  stack),  the  rate  of  change  in  9x/9x 
along  the  x  axis  is  negligible  and  the  term  3^x/9x^  is  zero.  A  further  sim¬ 
plification  of  the  diffusion  equation  is  obtained  if  a  continuously  emitting 
infinite  line  source  is  oriented  tiormal  to  the  mean  wind  direction.  Then, 
3x/3y  is  also  zero  and  Equation  4.16  is  reduced  to 


In  reality,  the  horizontal  wind  direction  is  rarely  steady,  source  emis¬ 
sions  vary  with  time,  and  line  sources  are  never  infinite.  When  using  a  model 
that  requires  these  assumptions,  it  is  prudent  to  have  information  on  the 
degree  of  departure  from  the  idealized  conditions  described  in  the  preceeding 
paragraphs.  The  information  required  involves  resolution  of  the  motion  fields 
that  the  plume  is  subjected  to  and  knowledge  of  the  plume's  source  geometry 
and  temporal  variability.  The  problem  is  further  complicated  by  the  fact  that 
although  the  plume  is  subjected  to  turbulence  as  it  is  carried  along  by  the 
wind  (Lagrangian  coordinate  system),  wind  and  turbulence  measurements  are  usu¬ 
ally  made  from  fixed  points  (Eulerian  coordinate  system).  Point  measurement 
devices  on  a  single  tower  are  rarely  adequate  for  providing  the  necessary 
field  measurements,  although  several  measurement  positions  along  the  path  of 
plume  travel  would  help  provide  the  required  information.  These  measurements, 
integrated  with  the  mesoscale  wind  field  measurements  in  a  well  defined  synop¬ 
tic  pattern,  are  required  to  obtain  an  adequate  description  of  atmospheric 
motion  fields. 

In  addition  to  difficulties  witti  K-theory  expressions,  modelers  are  faced 
with  the  problem  of  cloud  centroid  motions.  Labora'’  -y  measurements  (see,  for 
example,  Willis  and  Deardorff,  1976)  indicate  th-it  cloud  centroids  often  do 
not  remain  at  the  height  of  release.  Caughey  et  al .  (1983)  verified  this 
finding  in  their  convective  boundary  layer  (C3I.)  turbulence  structure  studies. 
Caughey  et  al .  note  that  downdrafts  cover  more  than  half  of  the  horizontal 
area  in  the  CBL,  the  effect  of  which  is  a  net  descent  of  material  released 
into  the  upper  part  of  the  CBL.  For  surface  releases,  the  centroid  can  either 
ascend  or  move  horizontally,  since  it  cannot  penetrate  the  ground.  In  convec¬ 
tion,  updrafts  lift  the  centroid  off  the  ground,  while  downdrafts  can  only 
cause  horizontal  spreading.  The  net  result  is  that  cloud  centroids  ascend 
from  the  surface  or  descend  from  aloft  until  they  reach  approximately  half  the 
boundary  layer  depth.  Consequently,  the  condition  w  =  0.  (Equation  4.15)  does 
not  guarantee  that  a  CBL  cloud  centroid  will  rei  ain  at  its  release  height. 
Measurements  of  divergence  field  changes,  vertical  wind  profiles,  and  atmo¬ 
spheric  stability  are  indicators  of  the  degree  of  convective  activity  and  the 
consequent  centroid  height  changes.  For  diffusion  tests,  downwind  sampling  to 
monitor  the  cloud  centroid  movements  is  nrudent.  Failure  to  identify  changes 
in  the  vertical  position  of  the  cloud  centroid  is  a  major  cause  of  uncertainty 
in  diffusion  cloud  characterization . 

Analytic  solutions  for  x  in  the  above  equations  follow  from  mathematical 
procedures  describing  molecular  conduction  of  heat  in  a  solid,  as  suggested  by 
Roberts  (1923).  Pasquill  and  Smith  (1983)  note  that  the  solutions  are  most 
readily  obtained  by  adopting  power  law  forms  of  the  K  and  wind  profiles, 

Kz(z)  =  Ki  (?^  )n,  u(z)  =  ui(^^)'n  (4.18) 


The  most  widely  used  analytic  solutions  for  diffusion  problems  are  the 
Gaussian  diffusion  equations.  These  are  obtained  using  the  following  assump- 
ti ons: 

a.  Exponents  m  =  n  =  0  in  Equation  4.18.  The  result  is  that  K  and  u  are 
taken  to  be  constant  with  height. 
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b.  Travel  distance  downwind  of  the  source  equals  ut,  where  u  is  a  trans¬ 
port  wind  and  t  is  travel  tine. 

c.  A  relationship  exists  between  K,  t,  and  the  standard  deviation  of 
plume  material  distribution  ( o)  so  that 


0  = 


(4.19) 


With  the  above  assumptions,  it  is  possible  to  take  advantage  of  the  fact  that 
the  area  under  the  Gaussian  (normal)  curve  is  unity.  Consider  the  z  component 
of  a  diffusing  cloud  moving  downwind  with  the  x  axis  defined  along  the  cloud 
centroid.  From  integral  tables. 


_J  F.XP  {-z^/2<4)  = 


(4.20) 


Then, 


_J  l/(/rir  az)  EXP(-zV2f^)  =  1.0 


(4.21) 


Similar  relations  hold  for  the  x  and  y  components.  Therefore,  the  distribu¬ 
tion  of  cloud  material  at  varying  distances  downwind  from  a  source  can  be 
described  as  a  function  of  cloud  sigmas,  providing  a  concise  statement  of 
continuity  for  use  in  diffusion  models. 

Given  tne  above  assumptions,  Gaussian  diffusion  equations  are  obtained  for 
the  appropriate  sets  of  boundary  conditions  (Pasquill  and  Smith,  1983): 

a.  Instantaneous  point  source.  The  initial  boundary  conditions  are: 


X  ■*■0  as  x,z 


(4.22) 


X  *  as  x,z  0. 


(4.23) 


K73x/9z  ■*■0.  as  z  ■►O.,  x>0. 


(4.24) 


I 


For  an  instantaneous  release  in  zero  wind,  the  source  strength  (Q  grams)  is 
the  total  volume  integral  of  concentration,  i.e., 


I  fix  AxAyAz  =  Q. 


(4.25) 


The  solution  for  instantaneous  point  source  Tickian  diffusion  from  Equation 
4.14  is; 


X  =  (Q/(2^(2TT)3/2axaya2))  EXP  [-  (^^/2ai  +  I2a\)] 


(4.26) 


b.  Continuous  point  source.  For  a  continuous  release  from  a  point 
source,  where  Q  is  expressed  in  units  of  grams/second,  the  additional  boundary 
conditions  required  to  obtain  a  solution  from  Equation  4.16  are: 


3x/3t  =  0..  32x/5x^  =  0 


(4.27) 


//xUAyAZ  -  q 


(4.28) 


The  solution  is; 


=  (Q/(27rava,  u))  FXF  [-  +  z2/2o|)] 


(4.29; 


c.  Continuous  infinite  crosswind  line  source.  The  additional  boundary 
condition  is: 


3x/3y  =  0. 


(4.30) 


For  the  continuous  release  from  a  line  source,  the  source  strength  Q 
(grams/meter  second)  is  equivalent  to  the  expression 


The  solution  to  an  infinite  line  source  obtained  from  Equation  4.17  is 


X=  [']n^  o^u))  EXP  [-  (z^n4)^ 


(4.32) 


The  three  solutions  given  above  do  not  include  provisions  for  calculation  of 
concentrations  at  points  not  at  the  cloud  centroid  or  for  reflection  of  the 
cloud  from  an  impervious  surface.  See  Sutton  (1953)  or  Pasquill  and  Smith 
(1933)  for  further  discussion  of  these  provisions. 

The  assjmptions  of  constant  u  and  K  are  very  restrictive  because  both  vary 
significantly  with  height.  To  compensate  for  thevariation  of  u  with  height, 
diffusion  models  make  use  of  a  "transport  wind",  u.  The  transport  wind  is  a 
mean  wind  obtained  by  averaging  over  the  depth  of  the  plume  and  for  plume 
travel  time.  This  wind  averaging  over  space  and  time  should  be  representative 
of  the  rate  of  the  diffusing  cloud's  downwind  travel.  A  number  of  procedures 
exist  for  estimating  u,  Cramer  et  ?!.  (1972)  of^er  an  algorithm  for  [i  based 
on  the  1  or  2  meter  wind  speed.  For  downwind  hazard  distance  calculations  at 
DPG,  the  first  1/2-  or  l-minute  elevationangle  readings  from  10-grams  ceiling 
balloons  are  frequently  used  to  estimate  u.  Continuous  time  and  space  aver¬ 
aged  vertical  wind  profiles  available  from  modern  remote  sensing  devices 
should  per:n't  more  satisfactory  u  estimates. 

To  compensate  for  assuming  a  constant  K  and  using  the  a=  /?Kt  relation, 
much  effort  is  expended  in  attempts  to  describe  o  as  a  function  of  atmospheric 
stability  parameters,  measured  turbulence,  and  downwind  t'^avel  distance.  The 
simplest  approaches  are  those  developed  by  Pasquill  (1961),  with  a  varying  as 
a  function  of  information  obtainable  from  general  weather  observations  (sun 
angle,  cloud  cover,  and  wind  speed).  The  Pasquill  procedure,  further  devel¬ 
oped  by  Turner  (1964)  is  currently  used  as  a  stability  category  input  for  DPG 
hazard  calculations.  An  updated  and  improved  form  of  this  procedure,  based  on 
net  radiation  measurement,  is  described  in  Myirski  (1983).  To  avoid  the  limi¬ 
tations  of  the  stability  category  approach,  measured  or  calculated  values  of 
turbulence  can  be  used  in  sigma  calculations  for  the  diffusion  equations. 

Sigma  e,  for  example,  is  derived  as  a  function  of  stability  parameter  z/L  in 
Section  7.  Sigma  z  can  then  be  calculated  as  Oq  x  f(x). 

K-theory  simplifications  leading  to  the  Gaussian  diffusion  equations 
permit  concise  solutions  for  diffusion  problems.  Given  a  source  strength, 
release  height,  configuration  and  duration,  cloud  dimensions  in  sigmas,  and  a 
transport  wind,  concentrations  downwind  from  a  source  can  be  readily  calcu¬ 
lated.  However,  the  performance  of  a  Gaussian  model  is  limited  and  expecta¬ 
tions  that  solutions  will  fit  a  given  set  of  measured  concentrations  cannot  be 
too  high  (Smith,  1984),  At  best,  many  replications  of  measurements  made  dur¬ 
ing  similar  atmospheric  conditions  will,  on  the  average,  approach  the  Gaussian 
calculated  value. 

Some  of  the  limitations  of  Gaussian  models  are  overcome  through  the  appli¬ 
cation  of  similarity  concepts  to  diffusion.  Similarity  models  have  been  based 
on  the  hypothesis  that  the  Eulerian  properties  which  characterize  the  SBL  can 
be  applied  to  the  Lagrangian  motions  of  particles  within  the  SBL.  Model 


levelopments  by  Horst  (1979)  and  Nieuwstadt  and  van  Ulden  (1978)  show  promise 
in  overcomini  some  Gaussian  model  limitations.  However,  similarity  theory  has 
not  successf  illy  treated  extremely  stable  nocturnal  situations  common  at  DPG. 
Even  with  this  limitation,  future  model  development  at  DPf  should  include 
consideratio  1  of  similarity  applications  to  diffusion.  New  K-theory  treat¬ 
ments  by  Wynjaard  and  Brost  (1984)  and  Fieoler  (1984)  should  also  be  investi¬ 
gated  for  applications  at  DPG. 

5.  KINETIC  ENERGY  AND  THE  RICHARDSON  NUMBERS 

Atmospheric  stability  is  related  to  turbulence  through  the  turbulent 
kinetic  energy  (TKE)  equation.  First  attempts  tc  relate  turbulence  to  stabil¬ 
ity  criteria  were  made  by  Mr.  L.  F.  Richardson  in  the  192()s  (Richardson, 

1925).  These  dimensionless  stability  criteria  came  to  be  known  as  Richardson 
Numbers.  A  complete  derivation  of  a  criterion  of  turbulence  for  the  atmo¬ 
sphere  was  made  by  Calder  (1949)  using  the  TKE  equation.  Cald^r  achieved  the 
derivation  by  defining  turbulent  kinetic  energy  a?  the  difference  between  the 
instantaneous  and  mean  kinetic  energy  equations,  with  particular  attention  to 
Reynolds  stress  terms.  Derivation  of  Richap^dson  numbers  from  the  TKE  equation 
in  this  section  follows  Calder's  procedure,  but  uses  simplified  notation. 

The  instantaneous  vector  equation  of  motion  is 


PdV/dt  +  pfc  =  ^0  -  W  +  uV^  (5.1) 

whei^e  dV/dt  is  a  total  derivative  of  the  vector  wind  while  f^-,  g,  Vp/p,  and 
pV^V/p  represent  forces  per  unit  mass  due  to  the  Coriolis  effect,  gravity, 
pressure  gradient,  and  viscoiis  dissipation.  The  dot  product  of  Equation  5.1 
with  the  velocity  vector  V  produces  the  kinetic  energy  equation  of  instantane¬ 
ous  motion 

'2pdV^/dt  =  -wgp  -  "V-VP  +  "V'pvV  (5.2) 


Note  that  the  Coriolis  term  drops  out  of  the  equation  because 

V -fc  =  X  V  =  0.  (5.3) 


Also,  because  the  force  of  gravity  only  operates  in  the  z  direction,  the  dot 
product  of  the  velocity  vector  with  this  term  reduces  to  the  scalar  wgp. 

On  the  left  side  of  Equation  5.2  the  dotted  velocity  vectors  have  been 
reduced  to  the  scalar  where 


=  (u^  +  +  w^). 


(5.4) 


The  total  differential  PdVVdt  expanded  into  its  component  parts  is 


pdV^/dt  =  p3V^/  3t  +  pV  -  W^ 


(5.5) 


In  the  next  series  of  steps,  the  total  differential  is  replaced  with  its 
components  in  '|he  desired  form.  The  continuity  equation  is  introduced  as  a 
product  with 


v^3p/3t  +  v^3(pu)3x  +  v^3(pv)/3y  +  V^3(pw)/3z  =  0.  (5.6) 


Because 


V^3p/  3t  =  3(V^p)  3t  -  p3V  Vst 


(5.7) 


and  using  similar  relationships  for  other  components  of  Equation  5.6,  the 
total  differential  in  Equation  5.2  can  be  replaced  by 


3(  Mf^)/ at  +  3(pV^u)/3x+  3((^/^)/^+  3((^^/)/32] 
=  -wg  p  -  V  ♦'5P  +  V 


(5.8) 


The  left  side  of  Equation  5.8  can  be  written  in  tensor  notation  following 
Calder  (1949),  yielding 

[3(PUi^)/3t+  3(  pui \)/ 3xk]  =  -  wgp  - 'V'VP  +  (5.9) 


The  instantaneous  equation  of  motion  is  next  separated  into  mean  ( p,  v) 
and  turbulent  (  p' ,  v')  components.  Because  the  average  of  each  turbulent 
component  is  zero,  the  average  value  of  terms  containing  a  single  turbulent 
component  is  zero.  However,  terms  containing  multiple  turbulent  components 
(p'u',  u'w',  etc.)  have  nonzero  averages.  Using  these  rules,  pu-j-  from  the 
left  side  of  Equation  5.9  becomes 

*wi  ^  =  ("p  +  "p')  (Hf  +  zITiiT^  +  TT^  ^)  =  \  ^  +  B'p'ITiU^  +  "p'u^^  (5.10) 
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The  terms  In  Equation  5.10  containing  p'  are  3  or  4  orders  of  magnitude 
smaller  than  terms  containing  p  and  are  usually  neglected.  The  second  term  on 
the  left  side  of  equation  5.9  is  treated  in  a  similar  manner.  The  result 
(neglecting  terms  with  p')  is 


pui^uk  =  p(ui^U|(  +  ui^U|(  +  2ui-ui‘u'|(  +  u^^u^)  (5.11) 


On  the  right  side  of  Equation  5.9,  wgp  is  expanded  to 


wgp  =  wgp  +  w'gp'  (5.12) 


In  this  case  w'gp'  is  not  neglected  because  it  will  be  the  remainder  after  the 
subtraction  of  two  larger  terms  in  Equation  5.19.  The  V*VP  term  is  expanded 
to 


^.7P  =  v^'VP  +  V'*VP'  (5.13) 

-i  o_\ 

The  final  term  on  the  right  side  of  Equation  5.9  is  V^pV'-V.  This  is  a 
dissipation  term.  With  expansion  into  mean  and  turbulent  components  it 
becomes  c  +  e',  where  e  represents  viscous  dissipation  into  heat. 

Gathering  terms  from  Equations  5.10  to  5.13  the  kinetic  energy  equation  of 
instantaneous  motion  is  given  by 


''2Pd(Ui"  +  u)^)/dt  =  i5[3(pUT^  +  pui^)/3t  +  3(p(Ui^U|(  +  U^^Uk  + 

_ _ ^  ^  __  (5.14) 

2u-ju'iu'|(  +  u'-j^ Uk ) )/3X|(]  =  -  wgp  -  w'gp'  -  V»VP  -  V'*VP'  +  e  +  e' 


The  vector  equation  of  mean  motion,  where  Equation  2.13  defines  the  u 
component,  is  used  as  the  starting  point  to  obtain  the  kinetic  energy  equation 
of  mean  motion.  The  major  difference  between  the  mean  and  instantaneous 
equations  is  the  treatment  of  Reynolds  stress  terms.  The  Reynolds  stress 
terms  are  represented  in  tensor  notation  as 


-pUjUk  =  TiK 


(5.15) 


Following  a  development  similar  to  the  instantaneous  equation  given  above,  the 
kinetic  energy  equation  of  mean  motion  is 
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'2oi(ui^)/dt  -  H[  3{  fUj  ^)/ 3t  +  3(  p(ui  ^k) )/ 3)(|(] 
=  -  Wrj  p  -  T  •'^  +  e  +  u-j  3T.jk/  3xk 


(5.16) 


Applying  the  rules  of  partial  differentiation  to  the  Reynolds  stress 
components  yields  the  expression 


Ui3Ti-k/3xk  =  3(TikUi)/3xk  -  Tik3ui/3xk 


(5.17) 


Reynolds  stress  components  are  also  found  in  the  kinetic  energy  equation  of 
instantaneous  motion.  From  Equation  5.14 


'2[  3(2  puiujui^)]/ axk  =  3(- TikUi )/ 3xk 


(5.18) 


Term  by  term  subtraction  of  the  kinetic  energy  equation  of  mean  motion  from 
Equation  5.14  for  instantaneous  motion  produces  Equation  5.19,  the  TKE 
equation. 


‘spcKlTi*  ^)/dt  =  3(  pu^'^)/ 3t  +  3(  pdTi  ^k))/ ^k3 

—  — 

=  -  wg  p  -  T'  ‘VP  +  e  +  Tjk  Sj,  /  3xk 


(5.19) 


The  term  w'gp'  is  transformed  into  more  useful  form  by  recognizing  from  the 
equation  of  state  that 


P7P  =  p7p  +  T7T. 


(5.20) 


Calder  (1949)  demonstrates  tiiat  the  ratio  of  pressure  fluctuations  to  mean 
pressure  is  several  orders  of  magnitude  smaller  than  corresponding  density  and 
temperature  ratios.  Also,  the  ratios  of  temperatures  and  potential  tempera¬ 
tures  (9)  are  nearly  equivalent.  The  approximation 


p7p  =  -  97  0 


(5.21) 


is  used  in  the  term  -w'gp'  to  produce 
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-w‘gp'  =  gp  w'e'/e  (5.22) 


The  remailing  terms  of  Equation  5.19  are  simplified  using  the  assumption 
of  horizontal  homogeneous  flow  along  the  direction  of  the  x  axis.  For  this 
condition,  mean  lateral  (v)  and  vertical  (w)  motions  are  zero  and  the  rate  of 
change  of  the  turbulent  components  in  the  horizontal  plane  is  zero.  Conse¬ 
quently,  the  term  ui^ u^  ,  where 

Uk  =  w  =  0  (5.23) 


is  zero  in  Equation  5.19.  The  term  ^3  (p  (u^^^ u^ )  )/3Xk  is  only  nonzero  for  the  -  ^ 

z  direction  is  shown  in  Equation  5.24  ^  J 


‘^3  (p  (u'-j^Uk  )  )/3Xk  =  ‘'23(p(u'^w'  +  v'^w'  +  w’^))/3z  (5.24) 


On  the  right  side  of  Equation  5.19,  V'*VP'  is  reduced  to  w'3P'/3z  because 
3P'/3x  and  3?'/3y  are  negligible  in  horizontally  homogeneous  conditions. 
Similarly,  Tikliui/3xk  becomes  -pu'w'3u/3z  because  v  and  w  are  zero. 

fathering  ter-ms  from  Equations  5.22  through  5.24,  the  mean  TKE  equation  for 
horizontal  homogeneous  iluw  is  thus 

[3(^i^)/3t]  =  ‘i[3(p(u'^  +  v'^  +  w’^))/3t  =  3(pE)/3t 

,  _ _  _  (5.25) 

-  -  +w'M)/3z  +  gp  w'e'/e  -  w'3P'/3z  -  pu'w'  3u/3  z  -  e ' , 

where  E  represents  turbulent  kinetic  energy.  If  a  further  assumption  is  made 
that  turbulent  conditions  are  stationary  (3E/3t  and  3p/3t  are  zero),  the  time- 
dependent  tet'ms  are  eliminated,  leaving  the  TKE  equation  composed  of  a  verti¬ 
cal  energy  flux  divergence  term,  a  buoyancy  energy  production  term,  a  term 
which  represents  energy  transfer  due  to  the  fluctuating  gradients  of  static 
pressure,  a  shear  energy  production  term,  and  viscous  dissipation.  The  par¬ 
tial  differentiation  rule  (as  applied  to  Reynolds  stress  components  in  Equa¬ 
tion  5.17)  is  applied  to  e'.  Viscous  dissipation  then  has  two  components, 

(u'^  v' )/"ix^  and  -  (v' 3^ u '  )/3 x^  and  comparable  terms  in  the  y  and  z 

directions.  Calder  (1949)  identifies  the  first  term  as  making  a  small  contri¬ 
bution  and  drops  the  term  from  the  equation.  The  second  term  represents  the 
mean  '"ate  of  dissipation  per  unit  volume  of  turbulent  energy  into  heat.  The 
sign  of  e'  in  Equation  5.25  was  changed  to  reflect  the  sign  of  the  second 
term. 
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The  flux  divergence  and  fluctuating  pressure  terms  are  not  amenable  to 
direct  measurement.  Consequently,  they  are  often  dropped  from  the  equations. 
Studies  by  Wyngaard  and  Cote  (1971)  and  Frenzen  (1983)  suggest  that  these 
terms  are  frequently  significant.  Frenzen  concludes  that  these  terms  con¬ 
tribute  to  the  TKE  equation  at  the  rate  of  one  fourth  of  the  sum  of  shear  and 
buoyancy  contributions  for  unstable  conditions.  Further  field  studies  of 
these  phenomena  are  required  to  clarify  the  role  of  these  terms  in  the  TKE 
budget. 

The  two  production  terms  in  the  TKE  equation  of  greatest  interest  are  the 
buoyancy  and  shear  production  terms.  The  buoyancy  production  term  can  be  pos¬ 
itive  or  negative  depending  on  the  sign  of  the  heat  flux  w'e'.  Vertical  heat 
flux  enters  the  TKE  equation  as  the  buoyant  force  per  unit  mass  (g/9)  that 
acts  on  an  eddy  whose  temperature  and  vertical  velocity  differ  from  the  sur¬ 
roundings  by  0'  and  w' .  As  the  eddy  travels  over  a  vertical  distance  i,  it 
performs  specific  work  proportional  to  (g/6)9't  that  enters  directly  into  the 
vertical  component  of  TKE.  On  the  other  hand,  the  shear  production  term  is 
always  positive  in  the  presence  of  the  stress  induced  by  air  movement  over  a 
rough  survade.  According  to  Calder  (1949),  it  represents  the  transformation 
of  mechanical  energy  present  in  the  mean  motion  to  the  mechanical  energy  of 
turbulent  motion.  The  shear  term  is  large  near  the  surface,  where  the  gradi¬ 
ents  3u/3z  are  large,  and  tends  to  decrease  with  height  in  the  SBL.  Shear 
usually  makes  a  larger  contribution  to  the  horizontal  than  the  vertical  compo¬ 
nent  of  TKE. 

Richardson  Numbers  are  formed  from  ratios  of  the  buoyancy  to  shear  produc¬ 
tion  terms  of  the  TKE  equation.  The  flux  Richardson  Number  is  a  dimensionless 
stability  parameter  defined  as 


Rf  =  (g/9’)(w'9' )/[(u'w’ )(3u/3z)]  (5.26) 


Rf  is  usually  not  calculated  from  direct  measurements  because  it  is  a  combi¬ 
nation  of  fluxes  and  gradients,  which  requires  substantial  amounts  of  instru¬ 
mentation.  A  more  commonly  used  form  of  Richardson  Number  is  the  gradient 
Richardson  Number  (Ri).  Ri  is  found  by  using  K-theory  flux-gradient  relation¬ 
ships  (Equations  4.8  and  4.10)  and  assuming  the  Kn,/Kh  ratio  to  be  unity. 

Ri  =  (Km/Kh)Rf  =  (g/0 )(90/3z)/(3u/3z)2  (5.27) 


_  Richardson  Numbers  are  zero  during  adiabatic  conditions  when  w'9‘  and 
39/9Z  approach  zero.  For  diabatic  conditions  with  an  upward  heat  flux  (due  to 
vertical  motion  and  lapse  of  potential  temperature  with  height  in  excess  of 
the  adiabatic  lapse  rate)  Richardson  Numbers  assume  negative  values.  The  con¬ 
verse  is  trie  for  a  downward  heat  flux  or  inversion  conditions. 

The  advantage  of  Ri  as  a  stability  parameter  is  that  the  gradients  (30 /3z) 
and  (3u/3z)  can  be  calculated  from  wind  and  temperature  profile  data,  but  gra¬ 
dient  calculations  also  have  disadvantages.  Gradients  of  wind  and  temperature 


are  small  in  relation  to  the  absolute  values  of  those  variables.  The  adia¬ 
batic  lapse  rate  is  0.01‘’C/m  and  measurements  of  that  order  of  accuracy  are 
necessary  for  accurate  gradient  calculations.  Most  temperature  instrument 
accuracies  do  not  exceed  0.1°C.  Comparable  accuracies  with  wind  measurements 
are  also  required.  The  problem  of  making  accurate  calculations  of  wi  from 
wind  gradient  measurements  is  further  complicated  because  the  squared  value  of 
(3u/8z)  in  the  denominator  of  Equation  5.27  magnifies  errors. 

Another  complication  with  the  use  of  Ricliardson  Numbers  arises  because  the 
gradients  of  wind  and  temperature  near  the  surface  are  height  dependent.  Gra¬ 
dients  computed  from  measurements  made  at  the  0.5  to  4  m  above  ground  level 
(AGL)  will  differ  considerably  from  gradients  computed  from  simultaneous 
measurements  made  at  2  and  8  m  AGL.  Therefore,  Richar-dson  Numbers  vary  with 
height,  and  the  data  from  which  Richardson  Numbers  are  derived  must  be  care¬ 
fully  specified  when  Richardson  Numbers  obtained  from  different  data  sets  are 
compared.  Also,  the  temperature  ano  wind  measurements  used  to  obtain  3u/3z 
and  30/3Z  in  Equation  5.27  must  be  made  over  identical  time  and  height  inter¬ 
vals.  Gradient  measurements  not  obtained  this  way  yield  unique  values  of 
Richardson  Numbers  which  are  useless  outside  the  dci.a  set  from  which  they  are 
computed.  The  Obukhov  length  stability  paraiueter,  described  in  Section  6, 
does  not  exhibit  a  dependence  on  measurement  height  and  can  be  used  in  place 
of  Ri. 

If  wind  measurements  are  available  at  only  one  level,  another  form  of 
Richardson  Number,  known  as  the  bulk  Richardson  Number  (BU),  can  be  used. 

BU  =  i00($/e){3?/3z)(z2/i?  )  (5.28) 


Nhile  dimensionally  correct,  bU  is  totally  empirical.  Because  the  temperature 
gradient  in  the  numerator  is  small  compared  to  u^ ,  BU  tends  to  be  small. 
Therefore,  the  calculated  value  of  BU  is  usually  multiplied  by  100  for  compar¬ 
ison  with  other  stability  parameters.  In  spite  of  its  empirical  origins,  BU 
performs  surprisingly  well  as  a  criterion  for  turbulence  conditions  in  the 
stable  atmosphere  where  Ri  fails  because  of  uncertainties  in  wind  gradients. 

Richardson  Numbers  are  useful  as  stability  parameters  when  the  assumptions 
made  during  simplification  of  the  TKE  equation  ere  not  seriously  violated. 

When  flow  is  not  horizontal  or  homogeneous,  or  when  significant  changes  in  the 
flow  occur  with  time,  discarded  terms  in  Uie  TKE  equation  become  significant 
and  simple  stability  parameter  esximates  fail  to  represent  existing  turbulence 
conditions.  Therefore,  stability  parameters  should  be  used  with  caution  as 
criteria  for  turbulence  during  periods  of  rapid  transition  such  as  sunrise  or 
sunset,  or  in  the  presence  of  transient  rnesoscale  systems.  The  usefulness  of 
stability  parameters  in  complex  terrain  or  near  surface  discontinuities  is 
also  open  to  question.  As  micrometeorologica )  measurement  capabilities 
improve,  new  stability  parameters  which  include  more  TKE  terms  will  be 
developed. 


6.  THE  OBoKHOV  LENGTH  AND  z/L. 


Obukhov  (1971)  recognized  the  requirement  for  a  more  fundamental  stability 
parameter  than  the  gradient  Richardson  number.  Assuming  that  conditions 
remain  homogeneous  in  space  and  stationary  in  time,  momentum  (Fni)  and  heat 
fluxes  (Ffi)  remain  nearly  constant  wthin  the  SBL. 


Fn,  =  p'u’w'  =  -tq  =  constant  (6.1) 


F|^  =  Cpp  w'9'  =  constant  (6.2) 


The  Cp  of  Equation  6.2  •"epresents  specific  heat  at  constant  pressure. 

Obukhov  observed  that  wind  shear  (;iu/3zj,  which  is  very  large  near  the  sur¬ 
face,  decreases  rapidly  with  height  and  reasoned  that  there  is  a  height  above 
the  ground  vjhere  the  dynamic  shear  stress  contributions  balance  the  buoyant 
contributions  to  turbulence.  Tliis  height  or  lerictn  scale  is  the  Obukhov 
length  (L).  Monin  and  Obukhov  (1954;  conclude  that  SBL  turbulence  can  be 
described  as  a  function  of  three  parameters;  (g/9),  u^^.,  and  w‘0‘.  The 
Obukhov  length  is  composed  of  a  combination  of  cliese  parameters  representing 
the  ratio  of  buoyancy  to  shear  terms  in  the  TKE  equation,  which  is  similar 
to  the  flux  Richardson  Number  (Equation  5.26).  With  the  aid  of  Equations  3.3 
and  3.11,  shear  is  described  in  terms  of  u*.  Then, 

.  _  (g/e)  (w'9'}  _  (c/o)(w'9')  _  z  k(g/e)(w'9')  _  z 

K'f  “  .jm — J-  —  —  >  —  .  I  .1  ■  *  ■  I  ■  ■  ■  ^  f  C  0\ 

(u'w' ) (3 u/3z)  -u*^  (u*i)m/kz)  <^^1  -u* 

The  sign  of  L  is  chosen  so  that  L  is  negative  for  unstable  thermal  stratifica¬ 
tion.  Hereafter  in  this  report,  the  Obukhov  length  proceeded  by  a  minus  sign 
(-L)  implies  heat  flux  directed  away  from  the  surface,  associated  with  a  lapse 
of  potential  temperature.  Because  (g/9),  w'6',  and  u,^  remain  nearly  con¬ 
stant  with  height  in  the  SBL,  the  magnitude  of  -L  changes  little  within  the 
first  few  Lens  of  meters  of  the  atmosphere.  This  is  in  contrast  to  the 
Richardson  'lumbers  that  can  vary  as  a  function  of  height.  The  Obukhov  length 
has  units  of  meters.  Therefore,  the  ratio  z/L  is  dimensionless.  Obukhov  in¬ 
tended  for  the  heigiit  above  the  surface  where  z/L  =  1  to  be  the  height  where 
balance  exists  between  the  buoyant  and  shear  contributions  to  turbulent  kinet¬ 
ic  energy,  and  defined  -L  as  the  height  of  the  sublayer  of  dynamic  turbu¬ 
lence.  This  balance  is  actually  ac'iieved  at  a  smaller  value  of  z/-L  (see 
Section  7). 

Because  it  depicts  the  relative  contributions  of  buoyancy  and  shear  to 
turbulence,  z/-L  functions  as  a  dimensionless  stability  parameter.  When  the 
buoyant  contribution  is  small,  the  magnitude  of  -L  is  large  with  respect  to 
the  height  of  the  SBL  and  z/-L  approaches  zero.  Turbulence  for  this  neutral 
condition  is  mechanical,  caused  by  the  effects  of  roughness  elements  and 
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surface  shear  stress  [In  the  early  moi'ninj  hours,  as  the  sun  is  rising,  the 
atmosphere  passes  briefly  through  this  condition].  As  the  ground  is  heated, 
w'o'  becomes  significant  and  buoyant  forces  ifeKe  a  contribution  to  turbu¬ 
lence.  Monin  and  Yaglom  (1971)  found  that  bu-.'-'ancy  effects  contriouting  to 
turbulence  as  zZ-L  exceeds  0.03.  With  increasing  w'G'  (or  as  z  increases  with 
w'9'  constant),  buoyancy  contributions  increase  until  the  point  is  reached 
where  the  contributions  of  buoyancy  and  shear  forces  to  turbulence  production 
are  in  balance.  This  occurs  as  z/-L  approaches  0.35.  As  z/-l  increases  be¬ 
yond  0.35,  buoyancy  becomes  progressively  mo''e  domi-'-ant  until  a  state  of  wind¬ 
less  convection  is  reached.  Windless  c.onvecti on,  characterized  by  intense 
surface  heating  with  light  and  variable  winds,  is  th*^  condition  whern?  shear- 
production  is  negligible  (u*  +  0.0).  W-indh'Ss  convection  is  frequently 
observed  during  the  middle  of  summer  d,ys  at  OPC. 

To  this  point,  the  contribution  of  latent  heat  effects  on  stability  have 
been  neglected.  Dyer  and  Hicks  (1^970)  accounved  for  latent  heat  effects  on  L 
as  in  Equation  6.4, 

-L  CppT  u*^/kg((  h  O.',-’  '  ])  (6.4} 


where  Fi-,  and  F-.  represent  fluxes  of  sensible  a.’.d  late.;t  heat.  At  its  max¬ 
imum,  F]  is  only  16  percent  of  F;^.  i'l  available  irom  humidity  flux  meas¬ 
urements,  the  F|  factor  should  be  considere.'.  At  DFG,  omission  of  F-; 
would  generally  result  in  a  negligible  e-rror. 

There  are  several  aavantages  to  using  Obiikhov  length  as  rs  stability  para¬ 
meter.  It  does  not  contain  a  mixture  of  flux  and  gradient  terms,  which  make 
Rf  difficult  to  calculate.  Kh/Km  is  not  assumed  constant  as  in  Ri .  Unlike 
Rf  and  Ri ,  L  is  independent  of  measurement  height.  The  main  assumptions  are 
the  same  as  those  required  for  Richardson  Munbe'S,  i.e.,  conditions  must  be 
reasonably  homogeneous  in  space  and  stationary  in  time.  The  energy  divergence 
and  pressure  terms  of  the  TKE  eguafion  are  stiil  ignored.  Equation  6.3  for  L 
has  the  disadvantage  that  it  is  compose'!  of  covariance  terms  (w'e’  and  u,*,) 
which  li-sjst,  be  deri-./ed  from  rreasurenent':  of  wind  and  temperature  fluctuations 
at  high  data  rates  using  sensiti/e  equipment.  Although  sonic  anemometer/ 
thermometers  are  capable  of  these  types  of  measurements,  z/-L  can  als;  .le 
exp*'v'.sed  in  a.  gradient  form  tha>-  pe^mirs  stability  computations  using  daca 
rom  i-.'ss  expensive  instruments.  A  gradient  forr.  of  zZ-L  is  obtained  through 
-se  of  Equations  4.10,  3.3,  u*  from  En'ration  3.1.  ,  and  defining  a  as  the 
ratio,  shown  by 

‘  '  'h/'m  -  -  { w '  0  '  )  /  (  i-  Z z  .' ]/ [  ( u '  w  ‘  )  /  r  9  u/ '9  z )  ] .  (6.5) 


ter  'jimpl  i  f  1  ca  ti  on  ,  the  result  is  Equacioi;  6.6. 


zZ  1,  =  -  igz  (  .': /9z )  1 -I  ■  V  z/z,.,) -  .  (tuio u/e z  ) ) 


(6.6) 
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9^-  tempera  :ure  scale  (degrees  Kelvin) 
u  -  coefficient  of  viscosity  (gram  centimeter"^  second"M 
p  -  density  (grams  centimeter"^) 

Oq,  ayy  -  vertical  and  horizontal  wind  angle  standard  deviation  (degrees  or 
radians) 

"  square  root  of  the  variance  in  horizontal  and  vertical  wind 
(meters  second"^  ) 

°x  "  standard  deviations  of  lateral,  vertical,  and  along  wind 
distributions  of  diffusing  material  (meters) 

Oxt  "  component  of  due  to  along-wind  fluctuations 

Tq  -  surface  stress  (gram  centimeter"^  second"^) 

T  -  stress  tensor  or  shear  stress  (gram  centimeter"^  second"^) 

‘I’m*  ‘J’h  "  limensionless  wind  shear  and  temperature  gradient. 

X  -  concentration  (milligrams  rneter"^ ) 

tp  -  diabatic  influence  function  (dimensionless) 
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Ri  -  gradient  Richardson  Number  (dimensionless) 

Ric  -  Criti:al  Richardson  Nunoer  (dimensionless) 

Rf  -  flux  Richardson  Number  (dimensionless) 
t  -  time  (seconds) 

T-  time-averaged  temperature  at  a  specific  height  (degrees  Kelvin) 

T’  -  fluctuating  temperature  component  (degrees  Kelvin) 

~  time-ay^raged  wind  speeds  along  the  x,  y,  and  z  axes  (meters 
second  ) 

u,  V,  w,  -  wind  speeds  along  the  x,y,z  axis  (meters  second  M 

u',  v' ,  w'  -  fluctuating  components  of  wind  speeds  along  the  x,  y,  and  z  axes 
(meters  second  ') 

u  -  transport  wind  speed  (meters  second"") 
u*  -  friction  velocity  (meters  second  ') 

V  -  vector  wind  (ui  +  vj  +  wk)  'meters  second  M 
X  -  along-wind  travel  distance  (meters) 
z  -  height  of  measurement  above  surface  (meters) 
zi  -  depth  of  the  convective  nixing  layer  (meters) 

Zq  -  roughness  length  (meters  or  centimeters) 
a  -  ratio  of  "p  to  K,;,  (dimensionless) 
c  -  turbulent  dissipation  rate  (meters"^  second  ^) 

fs.  ep  -  rates  of  turbulent  kinetic  energy-  addition  by  si  earing  and  buoyancy 
forces  (meters^  second  ’) 

Z  -  a  variable  relating  height  to  wind  speed 

A  -  delta  -  indicates  an  incremental  slice  or  section 

X  -  wavelength  (centimeters) 

'7-  del  operator  (i3/9x  +  j  3/ 9y  +  k  3/ Sz) 

h  -  potential  temperature  (degrees  Kelvin) 
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APPENDIX  A.  SYMBOLS 


A,  B  -  dime.'isionloss  constants  in  the  turbulent  kinetic  energy  equation 

BU  -  bulk  Richardson  Number  (dimensionless) 

Cn  -  refractive  index  structure  parameter  (meter" 

Cp  -  specific  heat  at  constant  pressure  (calories  gram  ^  degree 

D  -  dose  (milligrams) 

d  -  displacement  heigtit  (meters) 

dd  -  dosage  (milligram  -  minutes/m  ) 

E  -  turbulent  kinetic  energy  (meters  second  ) 

f(  )  -  functional  relationship  to  variable  inside  brackets 

fc  -  Coriolis  force  per  unit  mass  (meter  second  ) 

Ep.  F]  -  flux  of  sensible  and  latent  heat  (milliwatts  meter  ^) 

Fm  -  flux  of  momentum  (gram  centimeter  second"  ) 

Fx  -  viscous  force  per  unit  mass  (centimeter  second  ) 

i,  j,  k  -  unit  vectors  in  the  x,  y,  z  directions 

k  -  von  Karnan's  constant  (dimensionless) 

Kh>  ^''m»  *^'l  "  diffusivi  ty  ^for  heat, momentum,  and  moisture 
(meters  ^  second  ^) 

L  -  Obukhov  length  (meters) 

1  -  length  (meters) 

a  -  mixing  length  or  molecular  length  scale  (meters) 
n  -  exponent  on  power  law  wind  profile  (dimensionless) 

P  -  pressure  (millibars) 

?22  "  vertical  pressure  component  (millibars) 

0  -  source  strength  (variable  units) 

R  -  mixing  ratio  (milligrams/grams) 

RE  -  Reynolds  Number  (dimensionless) 


b.  Test  Data  Quality  Control.  The  similarity  theory  based  equations 
describing  og  as  a  function  of  gradients  obtained  from  profile  data  are 
quite  accurate.  With  high  quality  data  collected  at  a  site  of  known  roughness 
during  quasi-stationary  meteorological  conditions,  the  computed  turbulence 
should  correspond  closely  to  og  obtained  from  sonic  anemometer  measurements. 
Intercomparison  of  measurements  from  several  tower  sites  would  provide  a  check 
on  stationary  conditions.  Discrepancies  in  excess  of  20  percent  would  indi¬ 
cate  that  either  some  equipment  was  out  of  calibration  or  the  meteorological 
conditions  were  too  variable  for  collection  of  a  coherent  turbulence  data  set. 

c.  Posttest  Analysis:  The  aforementioned  instrumentation  provides  suffi¬ 
cient  data  for  a  fairly  complete  description  of  meteorological  conditions  in 
the  SBL.  The  description  begins  with  incoming/outgoing  long-  and  shortwave 
radiation  from  which  radiation  fluxes  can  be  determined.  Near  surface  scin¬ 
tillation  measurements  indicate  the  magnitude  of  small  scale  eddy  activity 
generated  by  radiative  fluxes  and  the  decay  of  larger  scale  eddies.  Eddy  ac¬ 
tivity  on  this  scale  deteriorates  the  performance  of  all  optical  instruments. 
Intermediate  turbulence  scales,  which  dominate  the  turbulence  spectrum  within 

a  few  meters  to  ten  meters  of  the  ground,  govern  the  dissipation  rate  of  sub¬ 
stances  released  into  the  SBL.  Spectrum  analysis  for  this  scale  of  turbu¬ 
lence,  along  with  momentum,  heat,  and  moisture  flux/gradient  calculations  pro¬ 
vide  the  detailed  information  required  for  understanding  cloud  growth  and  dis¬ 
sipation.  The  influence  of  larger  scale  eddies,  which  bodily  move  the  cloud 
without  contributing  significantly  to  its  dissipation,  is  obtained  from  calcu¬ 
lation  of  divergence,  analysis  of  acoustic  sounder  records,  and  measurement  of 
mesoscale  wind  fields.  These  measurements  must  then  be  described  within  the 
context  of  the  meso-  and  synoptic-scale  meteorological  events  taking  place. 

The  accumulation  of  these  data  should  provide  the  analyst  with  sufficient 
information  for  an  in-depth  analysis  of  meteorological  conditions  influencing 
a  test. 

d.  Model ing.  Because  L  is  constant  within  the  SBL,  profile  data  used  to 
calculate  L  provide  sufficient  information  for  model  computation  of  z/-L  over 
the  depth  of  the  SBL.  Because  oq  is  described  as  a  function  of  z/zq  and 
z/-L,  entire  turbulence  profiles  up  to  heights  of  20  or  30  meters  can  be 
achieved  using  16-meter  tower  profile  measurements  obtained  from  towers 
located  at  representative  downwind  positions.  Supplemented  by  deeper  acoustic 
sounder  wind  and  mixing  depth  profiles,  these  tower  measurements  could  produce 
the  detailed  information  required  for  more  accurate  diffusion  models.  To 
accomplish  this,  methodology  work  is  required  to  develop  better  relationships 
between  turbulence  calculations  or  measurements  and  the  diffusion  parameters 
02  and  Ox.  Further  methodology  studies  are  also  needed  to  test  alternatives 
to  Gaussian  diffusion  models. 
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SBL  similarity  arguments  apply  well  to  the  vertical  component  of  turbu¬ 
lence,  but  ere  not  as  useful  for  describing  the  horizontal  components.  Hori¬ 
zontal  wind  data  frequently  contain  a  relatively  large  amount  of  energy  at  low 
frequencies.  These  low  frequency  changes  scale  with  meteorological  phenomena 
occurring  above  the  SBL  (such  as  depth  of  the  convective  mixed  layer)  and  are 
also  subject  to  wind  direction  trends  and  terrain  influences.  No  well  estab¬ 
lished  algorithms  are  available  to  describe  these  effects.  Consequently, 
horizontal  turbulence  components  are  best  obtained  by  direct  measurement. 

Surface  layer  similarity  arguments  described  in  this  document  do  not  apply 
extremely  close  to  the  surface  or  outside  the  SBL.  The  turbulence  spectrum 
near  the  ground  is  distorted  by  the  proximity  of  a  fixed  boundary  which  breaks 
up  large  eddies.  This  effect  is  particularly  noticeable  for  low  frequency, 
high-amplitude  eddies  which  develop  during  unstable  thermal  stratification. 
Because  the  momentum  gradient  becomes  large  very  close  to  the  surface,  stabil¬ 
ity  near  th?  surface  approaches  neutral  regardless  of  the  intensity  of  surface 
heating.  Carratt  (1980)  has  developed  a  methodology  for  describing  near  sur¬ 
face  conditions,  but  no  conclusive  independent  testing  of  this  methodology  has 
yet  come  to  the  attention  of  the  author.  Considei-ao ie  emphasis  has  been 
placed  on  extending  similarity  concepts  ''eyona  the  SBL  through  the  depth  of 
the  planetary  boundary  layer.  More  work  is  also  required  in  this  area  before 
firm  solutions  will  be  available.  Future  instrumentation  needed  to  probe  the 
planetary  boundary  layer  and  beyond  will  consist  largely  of  remote  sensing 
instruments.  These  may  include  lidar  wind  finding  and  real-time  cloud  sam¬ 
pling  systems,  updated  algorithms  for  acoustic  sounder  temperature  profiles, 
and  FM-CW  radar  wind  sounding  systems.  Drones  may  also  be  used  to  obtain 
Lagrangian  (moving  with  the  diffusing  cloud)  as  opposed  to  Eulerian  (point 
measurement)  data. 

The  instruments  and  similarity  theory  arguments  in  this  document  can  pro¬ 
vide  a  vast  increase  over  present  OFG  turbulence  measurement  and  analysis 
capabilities.  With  appropriate  displays  and  data  reduction  programs,  this 
enhanced  capability  has  the  following  applications: 

a.  Test  Go/No-Go  Criteria.  Meteorologists  and  project  officers  can  use 
stability  information  as  a  quantitative  means  of  evaluating  the  suitability  of 
meteorological  conditions  for  diffusion  tests.  Gradients  of  wind  and  temper¬ 
ature  obtained  from  tower  profile  data  can  be  used  for  real-time  calculations 
of  z/-L,  to  assist  with  test  control.  For  tests  requiring  near  neutral  condi¬ 
tions,  a  z/-L  less  than  0.03  (-L  greater  than  67  for  z  =  2  m)  is  required. 

When  instability  up  to  free  convection  can  be  tolerated,  a  test  limitation  of 
0.35  for  z/-L  (-L  =  5.7  m)  can  be  specified.  Beyond  a  value  of  0.35  for  z/-L, 
buoyancy  forces  become  large  and  the  probability  of  erratic  cloud  behavior 
increases.  When  z/-L  exceeds  1.0  (-L  =  2  m)  buoyancy  forces  dominate  turbu¬ 
lence  as  mechanical  effeccs  steadily  decrease.  Further  increases  in  instabil¬ 
ity  create  a  condition  of  windless  convection  where  horizontal  and  vertical 
accelerations  are  large  in  comparison  to  the  mean  wind  field.  The  erratic 
looping  behavior  of  cloud  plumes  occurs  under  these  conditions.  Windless  con¬ 
vection  or  other  conditions  where  accelerations  become  large  in  comparison  to 
the  mean  wind  field  would  likely  be  unsuitable  for  diffusion  testing. 
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becomes  much  more  difficult  at  other  sites  where  significant  roughness  changes 
occur  near  the  test  area.  Achieving  temporally  representative  measurements  is 
more  difficult  because  of  the  short  duration  (on  the  order  of  one  to  several 
minutes)  of  some  DPG  smoke/obscurant  tests.  Representative  mean  values  of 
wind  and  temperature  can  be  obtained  for  such  short  times,  but  representative 
variances  and  covariances  require  measurements  for  about  10  minutes.  Triple 
products,  such  as  u'w'w'  found  in  the  TKE  equation,  would  required  even  longer 
averaging  times.  Wyngaard  et  al .  (1971)  and  Wyngaard  and  Clifford  (1978) 
offer  a  possible  solution  that  may  assist  in  iiiaking  temporally  representative 
measurements.  Scintillometers  provide  path-averaged  measurements  of  and 
wind  components.  Wyngaard  and  Clifford  suggest  that  statistically  stable 
samples  of  path-averaged  turbulence  data  can  be  collected  using  only  one  per¬ 
cent  of  the  averaging  time  needed  for  comparable  point  measurements.  A  meth¬ 
odology  could  be  developed  to  take  advantage  of  path  averaging  techniques  for 
meteorological  measurements  supporting  smoke/obscurant  testing. 

9.  CONCLUSIONS 

Useful  stability,  turbulence,  and  diffusion  nt as.'re.ments  can  be  nade  if 
the  measurements  are  made  at  a  reasonable  site,  with  appropriate  equipment  and 
proper  atmospheric  conditions.  A  reasonable  site  is  relatively  homogeneous, 
with  no  extraneous  obstacles  disturbing  the  turbulence  characteristics  of  the 
flow.  A  good  site  also  permits  an  accurate  determination  of  roughness  length, 
with  Zq  independent  of  wind  fetch  and  a  negligible  displacement  height  (d). 
Appropriate  equipment  includes  an  array  of  18-meter  tall  towers  instrumented 
for  wind  speed  and  direction  at  2,  4,  8,  ano  16  meters  A6L,  and  a  high- 
frequency  turbulence  instrument  and  Lyman-a  at  6  meters.  The  number  of  in¬ 
strumented  towers  required  varies  -'/ith  grid  size  and  the  degree  of  homoge¬ 
neity.  Except  under  extreme  diabatic  conditions,  each  tower  should  be  repre¬ 
sentative  of  ambient  conditions  within  a  several  hundred-meter  radius  at  a 
relatively  uniform  site  on  DPG.  Tower  measurements  supplemented  by  path- 
averaged  winds  for  divergence  and  C^  from  scintillometers,  vertical  wind  and 
mixing  layer  depths  from  an  acoustic  sounder,  and  radiation  measurements  from 
precision  spectral  pyranometers  and  pyrgeometers  would  provide  a  fairly  com¬ 
plete  set  of  turbulence  data,  assuming  tnat  the  data  were  collected  in  a  well 
defined  meso-  and  synoptic-scale  setting,  and  over  a  sufficient  time  interval 
to  be  representative  for  the  higher  moment  data  required.  Additional  informa¬ 
tion  on  microscale  pressure  variations  is  desirable,  but  difficult  to  obtain. 

The  similarity  theory  arguments  presented  ■’n  this  report  are  valid  only 
for  neutral  through  unstable  thermal  stratification.  A  corresponding  set  of 
equations  has  been  offered  for  the  stable  case,  but  these  equations  apply  only 
to  stabilities  less  than  the  critical  Richardson  Number  (Ric).  Oke  (1970) 
states  that  for  stability  beyond  Ric,  the  atmosphere  is  not  fully  turbulent 
and  there  may  be  no  definite  forms  for  profiles  of  wind  and  temperature.  With 
clear  skies  and  strong  radiation  losses,  nocturnal  stability  conditions  at 
PPG  almost  always  exceed  Ric.  Empirical  Formulations  sucn  as  the  bulk 
Richardson  Number  seem  to  work  as  well  as  anything  under  these  conditions. 
Consequently,  bulk  R’chardson  dumber  formulations  are  routinely  used  for 
nocturnal  diffusion  testing  at  DPG. 


of  the  amount  of  energy  entering  and  leaving  the  system.  These  measurements 
permit  calculation  of  the  radiation  balance.  Optical  instrument  tests  require 
these  radiation  measurements  and  they  are  also  useful  in  defining  stability 
categories  (Myirski,  1983). 

Measurements  on  the  meso  and  synoptic  scale  are  also  needed  to  describe 
the  setting  for  diffusion  experiments.  The  presence  of  fronts,  gravity  waves, 
and  related  phenomena  have  significant  influences  on  diffusing  clouds.  These 
scales  of  motion  are  responsible  for  the  transport  of  diffusing  material  and 
generate  lee-side  eddies  which  cause  diffusing  clouds  to  behave  in  unexpected 
ways.  Adequate  description  of  meso-  and  synoptic-scale  motions  requires  an 
expanded  meso-meteorological  network  and  a  satellite  image  receiver. 

Test  site  wind  and  temperature  gradients  used  for  computing  turbulence 
parameters  are  obtained  from  fixed  point  measurements  on  meteorological 
towers.  The  Project  Kansas  data  (Izumi,  1.971)  indicate  that  wind  and  tempera¬ 
ture  gradients  decrease  rapidly  with  height  above  the  surface.  Because  of  the 
rapidly  decreasing  gradients,  tower  profile  data  above  16  meters  contribute 
little  additional  information  on  the  conversion  of  mean  flow  to  turbulent 
kinetic  energy.  Consequently,  a  i6-mete'*  tower  instrumented  for  measuring 
wind  and  temperature  at  2,  4,  8,  and  16  meters  AGL  should  provide  sufficient 
information  for  gradient  computations. 

Gradients  are  computed  from  profile  data  averaged  over  several  tens  of 
minutes.  Consequently,  instrument  response  time  is  less  critical  to  good 
gradient  data  than  accuracy  and  repeatability.  Accuracies  approaching  0.01“C 
and  0,01  m/sec  are  desired  for  gradient  computations.  Such  accuracies  are 
approached  only  when  excellent  equipment  is  used  and  a  great  deal  of  care  is 
exercised  in  the  measurement  process.  Desirable  instrument  characteristics 
include  freedom  from  calibration  drift,  absence  of  hysteresis  effects,  and 
insensitivity  to  external  voltage  variations.  A  high-quality  thermal  shield 
and  aspiration  system  is  also  needed  for  the  temperature  measurements.  Free¬ 
dom  from  external  voltage  influences  is  best  achieved  with  sensors  that  trans¬ 
mit  data  as  frequencies  rather  than  voltages.  Well-designed  anemometer  chop¬ 
pers  and  quartz  crystal  thermometers  have  this  desirable  characteristic. 

Adequate  instrument  characteristics  are  necessary  but  not  sufficient  for 
obtaining  good  wind  and  temperature  gradient  data.  Wind  and  temperature  must 
be  measured  over  identical  time  and  height  intervals.  The  geometric  mean  of 
profile  measurement  heights  car  be  used  as  the  height  z  where  z/-L  is  calcu¬ 
lated.  The  2-meter  level  is  a  good  standard  level  for  the  lowest  profile  mea¬ 
surement  at  DPG  because  it  is  not  so  low  that  the  surface  influence  signifi¬ 
cantly  distorts  the  similarity  relationships.  Measurements  closer  to  the  sur¬ 
face  should  require  fast  response  anemometers.  Because  diabatic  influences 
usually  cause  gradients  to  be  nonlinear,  multiple  measurement  points,  spaced 
logarithmically  in  height,  are  desirable.  Thus,  good  gradient  data  could  be 
obtained  with  profile  measurements  2,  4,  8,  and  16  meters  AGL  with  a  geometric 
mean  z  of  5,66  meters. 

Having  the  above  mentioned  instruments  present  and  collecting  data  does 
not  guarantee  an  adequate  description  of  the  turbulence  field  for  diffusion 
tests.  The  measurements  must  also  be  temporally  and  spatially  representative. 
Achieving  spatial ly-representative  measurements  is  not  difficult  at  DPG,  be¬ 
cause  of  the  wide  expanses  of  nearly  uniform  terrain.  However,  the  instru¬ 
ments  must  not  be  shaded  by  man-made  obstacles.  Spatial  representativeness 


refractive  index  structure  parameter  (0^)  which  provides  d  measure  of  turbu¬ 
lent  activity  in  the  small  eddy  sizes  (scintillations).  Lawrence  et  al, 

(1970)  state  that  optical  scintillations  are  caused  primarily  by  irregulari¬ 
ties  in  the  refractive  index.  The  irregularities  are  variations  in  air  dens¬ 
ity  which  have  the  size  of  a  Fresnel  zone,  i.e.,  {X1V2  where  X  is  a  wavelength 
and  i  a  path  length.  For  visible  light  and  path  lengths  on  the  order  of  a 
kilometer,  eddy  size  scales  causing  scintillations  are  on  the  order  of  a 
centimeter. 

Details  of  eddy  activity  in  intermediate  ranges  can  be  obtained  by  a  sonic 
anemometer.  Unlike  the  LED  scintillometer,  the  sonic  anemometer  is  a  fixed- 
point  measurement  device.  A  typical  20-cm  path  length  between  probes  is  the 
limit  of  resolution  on  the  small  (high  frequency)  eddy  scale.  However,  in  the 
intermediate  size  ranges  which  most  strongly  influence  diffusion  on  the  scale 
of  interest  at  DPG,  the  sonic  anemometer  is  able  to  provide  detailed  informa¬ 
tion  on  the  fluctuating  wind  componants  (u',  v',  and  w' )  and  the  mean  vjind 
field,  uncompromised  by  the  mechanical  damping  or  overshoot  of  mechanical 
vanes.  This  detailed  measurement  of  the  turbulence  spectrum  is  riecessary  for 
analysis  of  the  energy  in  parts  of  the  turbulence  spectrum  of  interest  for 
smoke/obscurant  work.  Sonic  anemometers  provide  sufficiently  detailed  data 
for  analysis  of  the  variance  and  covariance  0''  wind  components.  This  is  par¬ 
ticularly  valuable  for  obtaining  the  w  components,  which  are  otherwise  diffi¬ 
cult  to  measure.  Difficulties  with  sonic  anemometers  include  high  cost,  the 
requirements  for  precision  orientation  and  leveling  to  prevent  cross-component 
contamination  and  their  electronic  complexity.  Although  discharge  anemometers 
may  perform  similar  measurements,  little  is  known  of  their  characteristics. 

To  obtain  moisture  flux,  a  Lyman-rx  humidity  system  should  be  operated  concur¬ 
rently  with  the  sonic  anemometer.  The  product  of  w'  with  the  high-speed 
humidity  measurements  will  provide  moisture  flux. 

The  Doppler  acoustic  echo  sounder  is  another  instrument  which  provides 
important  turbulence  information.  This  instrument  uses  the  Doppler  shift  of 
returning  sound  pulses  to  provide  mean  layer  wind  information  in  a  vertical 
profile  [and  the  square  root  of  the  variance  in  horizontal  and  vertical  wind 
velocity  (o^,  and  Oy^)].  Additional  information  can  also  be  obtained 
on  the  depth  of  stable  or  convective  boundary  layers  (zj)  the  height  of 
elevated  inversion  surfaces.  Convective  eddy  size  is  largely  a  function  of 
Zi.  Large-scale  horizontal  variations  in  wind  direction  which  contribute  to 
the  magnitude  of  data  also  scale  with  z-j  (Pasquill  and  Smith,  1983). 
Detailed  turbulence  (a^,  oy,  and  oyy)  information  is  best  obtained  over 
the  depth  of  a  few  hundred  meters  with  a  high-frequency  "mini"  acoustic 
sounder.  Mixing  height  and  elevated  inversion  surfaces  aloft  require  larger, 
more  powerful  low  frequency  units  reaching  heights  near  one  kilometer.  Other 
instrument  systems  useful  for  zi  measurements  are  described  in  Kaimal  et 
al.  (1982). 

For  a  more  complete  understanding  of  boundary  layer  turbulence,  it  is 
necessary  to  describe  the  meteorological  setting  in  which  turbulence  occurs. 
This  setting  includes  the  incoming/outgoing  energy  that  drives  this  turbu¬ 
lence,  which  is  the  long  and  shortwave  radiation  and  new  radiative  flux  as 
measured  by  a  set  of  pyranometers  and  pyrgeometers .  Operating  in  pairs,  with 
one  pyrgeometer  and  pyranometer  facing  up  and  the  other  facing  towards  the 
ground,  the  amount  of  long  and  shortwave  radiation  passing  through  a  hori¬ 
zontal  surface  is  unambiguously  measured.  This  provides  a  fundamental  measure 
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Algorichms  for  determining  along-wind  cloud  growth  (a^)  are  less  well 
developed  than  those  tor  Oy  and  Wilson  (1981)  states  that  except  for 

very  close  to  the  source  wnere  the  plume  is  small,  vertical  diffusion  and 
shear  advection  are  the  dominant  factors  determining  alonq-wind  dispersion. 
Wilson  offers  the  following  equation  for 


=  g^CO.OQ  (x  -r  ((h  +  0.5c^)lri((h  +  0.17  )/zo ) ) )  ^  +  ( '^t  ^  (7.10) 


where  h  is  the  effective  source  height  above  the  ground  in  meters  and  Oxt  is 
the  component  of  Ox  due  to  horizontal  wind  fluctuations  (u').  Further 
experimentation  is  required  to  refine  the  algorithms  for  (Equation  7.9) 
and  Ox  (Equation  7.10). 

8.  MEASUREMENT  REQUIREMENTS  AND  INSTRUMENTATION 

Adequate  description  of  turbulence  in  the  atmosphere  requires  knowledge  of 
energies  in  the  spectrum  of  eddy  sizes  ranging  from  the  height  of  the  upper 
boundary  (height  of  the  convective  mixed  layer  in  the  daytime)  to  the  scale  at 
which  viscous  dissipation  begins.  Measurements  of  incoming/outgoing  radiative 
fluxes  that  initiate  diabatic  effects  are  also  needed.  Although  no  single 
instrumented  tower  provides  all  of  these  measurements,  point  measurement 
devices  can  provide  a  great  deal  of  information  about  turbulence  fields. 

These  instruments  must  collect  data  over  adequate  time  and  space  averages  to 
provide  representative  measurements.  Wind  measurements  in  particular  should 
be  taken  from  multiple  locations  when  fixed  point  measurement  devices  such  as 
anemometers  and  vanes  are  used.  Path  averaging  devices  greatly  reduce  the 
instrument  deployment  requirements. 

An  impo'^tant  instrument  for  measurement  of  turbulent  structure  is  the  LED 
scintillometer  (Ochs  and  Ting-i-Wang,  1978).  Scintillometers  measure  path- 
averaged  wind  components  perpendicular  to  scintillometer  beam  paths.  A  path- 
averaged  horizontal  wind  field  can  be  resolved  with  two  beams  at  right  angles 
to  each  other.  Scintillometers  operate  over  path  lengths  from  300  to  1500 
meters,  ideal  for  most  diffusion  testing  at  DPG.  Scintillometer  path-averaged 
wind  components  measured  over  a  triangle  can  be  used  to  calculate  horizontal 
divergence  (Kjelass  and  Ochs,  1974).  For  an  equilateral  triangle,  the  equa¬ 
tion  becom.'S 


7^1  ‘V  =  1  (  Vq  +V  120  7  24o)/  ^  ~  3z» 


(8.1) 


wheri'  Vq.  V120.  V2ito  ^*"6  wind  components  perpendicular  to  triangle  sides 
orierited  east-v/est,  30  -210,  and  150  -330.  Length  of  the  triangle  side  is  1, 
and  A  the  enclosed  area.  Horizontal  divergence  is  related  to  accelerations  in 
the  vertical  wind  field.  Since  use  of  diffusion  equations  requires  the  tacit 
assumption  that  the  atmosphere  is  nondivergent,  divergence  measurements  may 
help  explain  discrepancies  between  calculations  and  measured  results  from  dif¬ 
fusion  experiments.  The  LED  scintillometer  also  provides  information  on  the 
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Table  7.2  presents  4ini,  B(z/-L),  and  ow/u*  for  representative  values  of 
z/-L. 


Table  7,2.  B(z/-L),  and  a„/u*  as  Functions  of  z/-L 

( from  Bi 1  toft,  1983) . 


B  z/-L) 


00 

1.30 

18 

1.29 

36 

1.33 

63 

1.41 

90 

1.48 

30 

1.72 

00 

2.74 

Equation  7.7,  used  in  a  ratio  with  the  logari tnmic  wind-profile  (Equation 
3.12),  satisfies  Equation  7.1  to  produce  a  method  for  computing  Og 

^  1.3  ((1  +  15  zZ-L)-^/*^  1.8  z/-L)^/^  =  ^  *  o 

u/u*  [In  (z/zq)  -^]/k  u 

Equation  7.8  describes  og  as  a  function  of  stability  parameter  z/-L,  height, 
and  roughness.  Following  Businger  (1973),  k  is  assigned  the  value  0.35.  This 
equation,  used  with  the  gradient  equation  for  z/-L  derived  earlier,  permits 
computation  of  Og  for  any  height  within  the  SBL  from  wind  speed  and  tempera¬ 
ture  profile  data. 

The  algorithm 


(7.9) 


is  less  well  established  than  the  one  for  Oy.  Biltoft  (1981)  used  Project 
Prairie  Grass  data  (Barad,  1958  a,b)  to  obtain  analytical  solutions  for 
with  f(x)  a  continuous  function  of  cloud  travel  distance  and  bulk  Richardson 
Number.  As  previously  mentioned  with  the  Cy  equation,  the  solutions  for 
a 2  apply  to  diffusion  over  10-minute  periods.  Adjustments  are  needed  for 
periods  which  differ  significantly  from  this.  Waldron  (1983)  has  successfully 
applied  these  analytical  solutions  on  DPG  diffusion  tests.  Alternative 
calculation  methods  are  presented  in  Pasquill  and  Smith  (1983).  Unfortu¬ 
nately,  Project  Prairie  Grass  data,  which  has  served  as  a  basis  for  most 
computation  methods,  did  not  include  adequate  Og  measurements  or  cloud  cen¬ 
troid  height  information  beyond  100  m.  Consequently,  values  computed  for 
these  data  are  speculative  beyond  100  m.  More  complete  data  sets  are  required 
for  a  better  definition  of  f(x). 


The  variance  in  vertical  wind  velocity  is  a  Reynolds  normal  stress  component 
of  the  TKl  equation  (Equation  5.25).  The  TKE  equation,  with  3(pE)/3t  assumed 
to  be  zero  and  neglecting  pressure  and  flux  divergence  terms,  is  reduced  to 


e ' /p  =  g  w'9'/6  -  u'w'3u/9z  "  ^  ^'b 


(7.2) 


Panofsky  and  McCormick  (1960)  used  this  version  of  the  TKE  equation  with 
e‘/p  defined  as  ay,^/zA^.  Solving  for  the  result  is  Equation  7.3. 


Ow  =  A[z(cs  B 


(7.3) 


A  and  B  are  dimensionless  constants  to  be  determined  below.  Constants  A  and 
B  account  for  the  fact  that  contributions  from  shearing  (eg)  and  buoyant 
(ep)  kinetic  energy  do  not  make  equally  efficient  contributions  to  the 
horizontal  and  vertical  components  turbulence.  With  the  aid  of  Equations 
3.3  and  3.11,  ej  is  described  in  Equation  7.4,  and  the  relationship  between 
e^/cp  and  flux  Richardson  Number  (from  Equation  6.3)  is  presented  in 
Equation  7.5. 

=  (xz)^  {du/dz)^/4-i^^  (7.4) 


Equations  7.4  and  7.5  are  then  applied  to  Equation  7.3.  After  some  algebraic 
manipulation  and  division  by  u*  the  result  is  Equation  7.6 


Ovy/u*  =  Ak”’ ((tm  +  B(z/-L))^/^ 


(7.6) 


Following  the  findings  of  Biltoft  (1983),  constants  Ak  and  B  were 
assigned  values  of  1.3  and  1.8,  respectively. 


The  dimensionless  wind  shear  term  ((fpi)  in  Equation  7.6  is  a  slowly 
decreasing  function  of  z/-L,  whereas  the  buoyancy  term  (B  z/-L)  is  an 
increasing  function  of  z/-L  (Table  7.2).  The  combined  effect  of  these  terms 
is  that  oy,/u*  remains  nearly  constant  at  1.3  until  z/-L  reaches  a  value 
of  0.35.  This  is  the  equilibrium  point  w*-  shear  (ijim)  is  in  balance  with 
buoyancy  (1.8  z/-L).  Beyond  this  point,  buoyancy  becomes  dominant  and 
increases  rapidly  with  z/-L.  Using  Equation  6.8  for  the 
working  expression  for  Equation  7.6  is 


o„/u*  =  1.3((1  +  15  zZ-L)-'/**  +  1.8  z/-L)^/^  (7.7) 
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described  by  Irwin  (1979)  as  a  complex  function  of  time,  Lagrangian/Eulerian 
time  scales,  stability,  mixing  layer  depth,  release  height,  and  friction 
veloci ty . 


Table  7.1,  Function  f{x)  for  Indicateo  Downwind  Travel 


Di stances 

and 

Stabi  1  i  ty 

Condi tions. 

Downwi nd 

Travel  Distance  From  Point  of 
(x,  iiioters) 

Release 

Stabi  1  i  ty 

100 

200 

400 

1000  2000 

3000 

Unstable 

0.95 

0.87 

0.7 

0.48  0.37 

0.3 

Neutral 

0.8 

0.7 

0.65 

0.42  0.35 

0.3 

Stable 

0.64 

0.51 

0.45 

0  3'  0.32 

0.3 

Sigma  A  measurements  required  for  use  with  the  factors  in  Table  7.1  should  be 
obtained  from  data  collected  at  a  rate  of  one  per  second  over  10  minutes.  For 
measurements  over  times  significantly  different  from  10  minutes,  adjustments 
are  required  (see,  for  example,  Cramer  et  al . ,  1972). 

Attempts  to  apply  similarity  arguments  to  horizontal  wind  components  have 
not  worked  well.  Horizontal  wind  angle  data  are  subject  to  mesoscale  circula¬ 
tions  and  local  terrain  influences  which  impose  trends  on  the  data.  Direct 
measurements  of  are  best.  Fortunately,  on  varies  little  with  height 
within  the  constant  shear  stress  layer  and  the  required  horizontal  wind  angle 
measurements  are  readily  obtained  from  high  quality  micrometeorol ogical  wind 
vanes . 

The  Gaussian  diffusion  equation  describes  the  vertical  growth  of  a  cloud 
ia^)  as  a  function  of  the  vertical  turbulence  component.  This  component  is 
often  represented  as  og,  the  standard  deviation  of  vertical  wind  angle  fluc¬ 
tuations  (frequently  obtained  from  biaireccional  vane  measurements).  Experi¬ 
ence  at  DPG  with  bidirectional  vane  measurements  indicates  that  these  instru¬ 
ments  usually  operate  well  when  positioned  A  or  more  meters  above  the  surface, 
but  sometimes  fail  to  produce  representative  data  under  extreme  diabatic 
condi tions. 

In  addition  to  direct  measurement,  it  is  possible  to  calculate  og  as  a 
function  of  stability  using  similarity  th?,ory  i^elationships.  Sigma  e  is 
approximated  by  the  ratio  of  (the  square  root  of  the  variance  in  vertical 
wind  velocity,  w'w' )  to  u,  as  in  Equation  7.1. 


Og  =  o^/ u  (7.1) 
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Equatloi  6.6  offers  several  advantages  as  a  means  of  describing  stability 
in  the  SBL.  First,  it  does  not  contain  k.  The  absolute  value  of  von  Karman's 
constant  is  not  precisely  known  isee  Businger,  1973);  accepted  values  range 
from  0.34  to  0.41.  Second,  Equation  6.6  contains  no  high  order  gradient  terms 
such  as  (3ii  3z)^  in  Equation  5.27  for  Ri .  This  is  an  important  consideration 
because  the  gradient  terms  are  the  largest  source  of  uncertainty  in  Ri 
expressions . 

Equation  6.6  requires  a  specification  of  a  for  solution.  The  term  a, 
shown  in  Equation  6.5  as  the  ratio  Kf,/Kni,  is  identified  by  Businger  et  al . 
(1971)  as  the  ratio  of  the  dimensionless  wind  shear  (Equation  3.11)  to  the 
dimensionless  temperature  gradient  (-tj,) 


‘  =  K^/Km  =  'tm/'J’h  =  C(kz/u*)/(3u/9z)]/[(kz/9*)/(3e/82)]  (6.7) 


where  0^^  is  a  temperature  scale  analogous  to  u*.  Businger  et  al.  (1971) 
used  Project  Kansas  data  (Izumi,  1971)  to  establish  empirical  relationships 
between  ijifp,  1>h>  ^ihd  z/-L,  as  represented  in  Equations  6.8  and  6.9 

=  (1  15  z/-!.)-^/**  (6.8) 


<^h  =  0.74(1+  9  z/-L)'^/^  (6.9) 


Based  on  Equations  6.7,  6.8,  and  6.9,  a  is  specified  as  a  function  of  z/-L  by 
the  expression 

t  =  1.35(1  +9  z/-L)^/2/(i  +  15  z/-L)l/-+  =  (5,10) 


With  a  and  (Equation  3.13)  specified  as  functions  of  z/-L,  a  simple  conver¬ 
gence  scheni!  can  generate  rapid  computer  solutions  for  z/-L  by  using  wind  and 
temperature  gradient  data. 

7.  STABILITY,  TURBULENCE,  AND  DIFFUSION 

As  stat'd  in  the  introduction  and  further  described  in  Section  4,  Gaussian 
'iittjsi  .n  •  , nations  rinjuire  specification  of  cloud  growth  sigmas,  Pasquill 
(!97t),  *i‘  f  arther  oiabo’^ation  by  Doran  et  al.  (1978),  used  field  measure- 
'Tnnt.  *  ,  »  ’/  y  as  a  tunctior  of  horizontal  turbulence  components  (o/\), 

trav  ’  ■  stm  >,  and  a  distance  dependent  factor  f(x)  as  in  Equation  1.2. 
^nn^tiy  presented  in  Table  7.1,  based  on  the  findings  of  Pasquill 

(It/e  V  :  ri  .-t  al.  (19781.  At  distances  beyond  1000  m,  the  neutral  case 
fix;  vi’e  «.r.  taken  to  be  intermediate  between  the  unstable  and  stable 
valu'  , ,  asf  pt)tically  approaching  0.3  at  extended  distances.  These  numbers 
represent  ,  imp! i tied  estimates  of  turbulence  effects.  Function  F(x)  is 
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APPENDIX  C.  ABBREVIATIONS 


AGL  -  above  ground  level 

CBL  -  convective  boundary  layer 

DPG  -  US  Army  Dugway  Proving  Ground 

NIOSH  -  National  Institute  of  Occupational  Safety  and  Health 
OSHA  -  Occupational  Safety  and  Health  Act 
SBL  -  surface  boundary  layer 
TKE  -  turbulent  kinetic  energy 
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